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ABSTRACT 

We analyse the observed fractions of core-collapse supernova (SN) types from the 
Lick Observatory SN Search (LOSS), and we discuss the corresponding implications 
for massive star evolution. For a standard initial mass function, observed fractions 
of SN types cannot be reconciled with the expectations of single-star evolution. The 
mass range of Wolf-Rayet (WR) stars that shed their hydrogen envelopes via their own 
mass loss accounts for less than half of the observed fraction of Type Ibc supernovae 
(SNe Ibc) . The true progenitors of SNe Ibc must extend to a much lower range of initial 
masses than classical WR stars, and we argue that most SN Ibc and SN lib progenitors 
must arise from binary Roche-lobe overflow. In this scenario, SNe Ic would still trace 
higher initial mass and metallicity, because line-driven winds in the WR stage remove 
the helium layer and propel the transition from SN lb to Ic. Less massive progenitors of 
SNe lb and lib may not be classical WR stars; they may be underluminous with weak 
winds, possibly hidden by overluminous mass-gainer companions that could appear 
as B[e] supergiants or related objects having aspherical circumstellar material. The 
remaining SN types (II-P, II-L, and Iln) need to be redistributed across the full range 
of initial masses, so that even some very massive single stars retain H envelopes until 
explosion. We consider the possibility of direct collapse to black holes without visible 
SNe, but find this hypothesis difficult to accommodate in most scenarios. Major areas 
of remaining uncertainty are (1) the detailed influence of binary separation, rotation, 
and metallicity, (2) mass differences in progenitors of SNe Iln compared to SNe II- 
L and II-P, and (3) the fraction of SNe Ic arising from single stars with the help 
of eruptive mass loss, how this depends on metallicity, and how it relates to diversity 
within the SN Ic subclass. Continued studies of progenitor stars and their environments 
in nearby galaxies, accounting for SN types, may eventually test these ideas. 

Key words: binaries: general — stars: evolution — stars: mass loss — supernovae: 
general 



1 INTRODUCTION 

The observed fractions of various types of core-collapse su- 
pernovae (CCSNe) provide key information about the evo- 
lution and ultimate fates of massive stars. Because of their 
tremendous luminosity, SNe can potentially be used as di- 
agnostics of mass loss and the evolution of individual stars 
at great distances and in a variety of galactic environments, 
but only if we first understand how to map initial masses and 
evolution of different progenitor stars to the various types of 
SNe that they produce. If drawn from a stellar population 
that obeys a standard initial mass function (IMF), the ob- 
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served fractions of different CCSN types constrain the ranges 
of initial mass for their progenitors, as well as the evolution- 
ary paths they take before death. The aim of this paper 
is to explore how the IMF can be sampled in order to be 
consistent with the observed fractions of CCSN types. 

The main observed types of SNe that we consider are 
II-P (plateau), II-L (linear), Iln (relatively narrow lines), 
lib (transitional), lb, Ic, and Ibc-pec (see §2). Spectroscopic 
classification criteria for these are reviewed by Filippenko 
(1997). Pre-SN mass loss of the progenitor star determines 
which of these types of SN is seen, stripping away various 
amounts of the star's H and possibly He envelopes before 
core collapse ejects the remaining envelope. SNe Ibc are the 
result of complete removal of the H envelope, SNe lib have 
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retained only a small H mass (typically <0.5 M©), while 
SNe II-L and II-P have retained increasingly more of their 
H envelopes. SNe Iln are different in the sense that their 
spectral appearance is determined largely by shock interac- 
tion with circumstcUar material (CSM) lost in the decades 
or centuries preceding core collapse. 

The three potential mechanisms for a SN progenitor's 
mass loss are via steady winds, eruptive mass loss, or mass 
transfer due to Roche-lobe overflow (RLOF) in a close bi- 
nary system. Depending on which dominates, the amount of 
mass lost may or may not depend in a simple way on metal- 
licity or on the initial mass of the progenitor star, making 
reliable predictions diflicult without a more complete under- 
standing of mass loss. The evolutionary state — red super- 
giant (RSG), blue suporgiant (BSG), luminous blue variable 
(LBV), and Wolf-Rayet (WR) stars of the WN and WC se- 
quences — and hence the stellar radius at the time of explo- 
sion, are also important, although these can be considered as 
largely the result of mass loss. Massive stars have substantial 
steady stellar winds through most of their lives (see Lamers 
& Cassinelli 1999), with either metallicity-dependent, line- 
driven winds in hot stars (Kudritzki & Puis 2000), or slow, 
pulsation/dust-driven winds in cool stars (Rcimcrs 1977). 
Single-star evolution models adopt simple prescriptions for 
these steady winds (Meynet et al. 1994), and aim to predict 
the fates of massive stars as functions of initial mass and 
metallicity (e.g., Heger et al. 2003; Eldridge & Tout 2004). 

Recent observational work, however, has demonstrated 
that the standard observational mass-loss rates used as input 
to these models are far too high; the standard mass-loss 
rates of hot stars (e.g., Nieuwenhuijzen & de Jager 1990; de 
Jager et al. 1998) arc reduced by factors of 3-10 when the 
effects of clumping are considered properly (Bouret et al. 
2005; FuUerton et al. 2006; Puis et al. 2006). There is also 
a parallel problem in cool star mass-loss rates — reduction 
of an order of magnitude to the standard Reimers formula 
for rod giants may be required (Meszaros et al. 2009), and 
it would be interesting if this also affects more massive RSG 
stars. These lower mass-loss rates have a profound impact 
on stellar evolution and SN progenitors, requiring us to turn 
to cither eruptive mass loss (Smith & Owocki 2006) or close 
binaries (e.g., Paczyriski 1967) to make up the deficit. We 
will see that this turns out to be a major theme in explaining 
the frequencies of SN types. 

The stripping of a star's H envelope due to mass trans- 
fer in RLOF binary systems has long been considered a 
likely mechanism to produce WR stars and the progeni- 
tors of SNe Ibc (e.g., Paczynski 1967; Podsiadlowski et al. 
1992). Recent stellar evolution models attempt to account 
for this (e.g., Eldridge et al. 2008), but considerable uncer- 
tainty surrounds empirical estimates of binary fractions; see 
Kobulnicky & Fryer (2007) and references therein. 

In addition to close binary evolution, a major uncer- 
tainty concerns the net effect of episodic and eruptive mass 
loss during late stages of stellar evolution (Smith & Owocki 
2006). These outbursts are observed to occur, and studies 
suggest that they shed more mass from a star than do steady 
winds (Smith & Owocki 2006). The importance of sudden, 
short-duration eruptive mass loss is a concern for the pre- 
dictive power of any stellar evolution model, none of which 
currently include it. Observational clues from CSM inter- 
action in SNe Iln dictate that heavy mass loss sometimes 



occurs shortly before core collapse (e.g., Smith et al. 2007, 
2008b, 2010; Chugai et al. 2004); if heavy mass loss is con- 
centrated in brief events during the last few thousand years 
before core collapse, then the statistical distribution of end 
fates (i.e., SN types) won't necessarily reflect the observed 
relative fractions of WN, WC, RSG, LBVs, and so on, which 
are determined by the time spent in each state. This is crit- 
ical and potentially misleading, since many stellar evolution 
codes are linked to these observed fractions. 

Another key point is that both binaries and eruptions 
are probably less sensitive to metallicity than line-driven 
winds of hot stars. Some studies have shown that the ob- 
served fraction of SNc Ibc compared to SNe II increases 
with metallicity, implying that metallicity-dependent winds 
play an important role (Prantzos & Boissier 2003; Prieto 
et al. 2008; Boissier & Prantzos 2009). On the other hand, 
observations have also revealed a large population of WR 
stars in low-nictallicity galaxies, which cannot be explained 
by stellar winds alone (Izotov ct al. 1997; Brown et al. 2002; 
Crowther & Hadfield 2006), while the broad- lined SNe Ic 
that accompany gamma-ray bursts (see Woosley & Bloom 
2006 for a review) seem to prefer low metallicity (Stanek et 
al. 2006; Modjaz et al. 2008). 

There have been a few previous investigations of rel- 
ative SN rates that our study builds upon. Cappellaro et 
al. (1997) examined the statistics from 110 SNe (including 
SNe la), deriving widely adopted rates of various SN types in 
different environments. More recently, Smartt ct al. (2009) 
considered a volume-limited sample of nearby CCSNe and 
examined the relative fractions of SNe II-P, II-L, lib, Iln, 
lb, and Ic, as we do. Our results arc different from theirs, 
as described below, leading to some quite different impli- 
cations for massive stars. Finally, Arcavi et al. (2010) have 
recently submitted a paper independent of our study using 
SNe from the Palomar Transient Factory (PTF), finding a 
difference in the relative fractions of SNe II, lib, lb, Ic, and 
broad-lined Ic between large galaxies and dwarf galaxies. 
Since our survey did not adequately sample dwarf galaxies, 
the study by Arcavi et al. (2010) is complementary to ours, 
although we find significantly different relative fractions of 
SN types in large galaxies. We also consider direct detections 
of SN progenitor stars from pre-explosion data, discussed in 
considerable detail below and reviewed recently by Smartt 
(2009). Throughout, we include this information along with 
current ideas about massive single and binary stars. 

Here we present and discuss the implications of the ob- 
served relative fractions of different types of CCSNe in a new 
volume-limited sample, measured during the course of the 
Lick Observatory Supernova Search (LOSS) conducted with 
the Katzman Automatic Imaging Telescope (KAIT; Filip- 
penko et al. 2001). This follows a series of papers discussing 
LOSS. Paper I (Leaman et al. 2010) describes the method of 
deriving rates from LOSS data. Paper II (Li et al. 2010a) dis- 
cusses the luminosity functions of SNc and gives a detailed 
discussion of how the different fractions of SN types were 
derived, and Paper III (Li et al. 2010b) presents relations 
with host galaxies and other details. 
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Table 1. Volume-limited core-collapse SN fractions 



SN Type 


fraction 


error 




( % ) 


{ % ) 


Ic 


14.9 


+4.2/-3.8 


lb 


7.1 


+3.1/-2.6 


Ibc-pcc 


4.0 


+2.0/-2.4 


lib 


10.6 


+3.6/-3.1 


Iln 


8.8 


+3.3/-2.9 


II-L 


6.4 


+2.9/-2.5 


II-P 


48.2 


+5.7/-5.6 


Ibc (all) 


26.0 


+5.1/-4.8 


Ibc+IIb 


36.5 


+5.5/-5.4 
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Figure 1. Relative fractions of CCSN types in a volume-limited 
sample from LOSS. This is slightly different from the fractions 
quoted in Paper II, in order to better suit the aim of this paper 
as explained in the text. The main difference is that we exclude 
SNe in highly inclined galaxies because of extinction effects, and 
we reorganise the class of SNe Ibc-pec (namely, we moved broad- 
lined SNe Ic from the "Ibc-pec" category to the "Ic" group). 



2 OBSERVED CCSN FRACTIONS 

Figure [T] shows a pie chart illustrating the relative fractions 
of different types of CCSNe derived from LOSS. These val- 
ues are taken from the volume-limited fractions of all SN 
types derived in Paper II, with the thermonuclear (Type la) 
explosions subtracted from the sample. The relative frac- 
tions of the total for CCSNe are listed in Table 1, and these 
values are adopted throughout this work. See Paper II for 
further details on how these numbers are derived from our 
survey. Errors in Table 1 were estimated using a random 
Poisson number generator to sample from a list of fake SNe 
with fractions corrected for various observing biases, with 
lO'' realizations. Paper II discusses this in more detail. 

There are several important points to note here. This 
volume-limited sample of CCSNe excludes most of the 
so-called "SN impostors" (e.g.. Van Dyk 2010; Smith et 



al. 2010, in preparation), which appear as relatively faint 
SNe Iln that are often discovered by KAIT. If we had in- 
cluded them, the fraction of SNe Iln would be significantly 
higher; note that even without the SN impostors, however, 
our relative fraction of SNe Iln is higher than in previous 
studies (CappeUaro et al. 1999; Smartt 2009). The crite- 
ria for excluding an individual SN impostor are admittedly 
somewhat subjective, but this is a necessary step since the 
diversity and potential overlap of SNe Iln and massive star 
eruptions are not fully understood yet. Generally, if an ob- 
ject has a peak absolute R or unfiltered magnitude brighter 
than —15 and has line widths indicating expansion speeds 
faster than about 1000 km s~^, we include it as a real SN Iln. 
Less luminous and slower objects are considered impostors 
and are excluded. 

Unlike previous studies, we include a category called 
"SNe Ibc-pec" (peculiar; see Paper II). This category was 
necessary to introduce in Paper II because some SN Ibc 
vary significantly from the template light curves used to de- 
rive the control times for SNe lb and Ic. As such, the "Ibc- 
pec" category in Paper II includes some broad-lined SNe Ic 
such as SN 2002ap that are clearly SNe Ic. We have moved 
these to the SN Ic category for the purpose of this paper, 
since they clearly correspond to massive stars that have fully 
shed their H and He envelopes. This has a small effect on the 
overall statistics, because broad-lined SNe Ic are very rare in 
our sample, contributing only 1-2% of all CCSNe. This is in 
agreement with the recent study of Arcavi et al. (2010), who 
find that broad-hned SNe Ic contribute only 1.8% of CCSNe 
in large galaxies. It is noteworthy, however, that Arcavi et 
al. (2010) find broad-lined SNe Ic to be much more common 
(~13% of CCSNe) in low-metallicity dwarf host galaxies. 
We also exclude SNe occurring in highly inclined galaxies, 
where dust obscuration may introduce statistical problems 
that are difficult to correct. As a result of these minor adjust- 
ments, made because our goal of investigating implications 
for massive-star evolution is different from the goal of deriv- 
ing relative rates and correcting for observational biases, the 
relative fractions of various SN types in Table 1 and Figure[T] 
differ slightly from the results in Paper II. 

In quoting fractions of various SN types, we ignore 
metallicity, galaxy class, and other properties, although we 
are cognizant of the importance of these properties and con- 
sider them in our discussion below. The galaxies included in 
the LOSS survey span a range of luminosity, with most of the 
CCSN hosts corresponding roughly to metallicities of 0.5-2 
Zq (Garnett 2002; the LOSS galaxy sample spans a range 
of Mk from about —20 to —26 mag, but most of the CCSN 
hosts are in the range —22 to —25 mag; see Paper II). We 
note some trends in Paper II, such as the fact that SNe Iln 
appear to prefer lower luminosity spirals, whereas SNe Ibc 
seem to prefer large galaxies and therefore higher metallicity, 
consistent with previous studies (Prantzos & Boissier 2003; 
Prieto et al. 2008; Boissier & Prantzos 2009). LOSS is biased 
against very faint dwarf galaxies, since larger galaxies with 
potentially more SNe were targeted to yield a richer harvest 
of SNe. However, low-luminosity galaxies seem to have more 
than their expected share of star formation per unit mass, 
and probably contribute 5-20% of the local star formation 
(Young et al. 2008). If unusually luminous SNe Iln and II-L 
favour such low-luminosity galaxies, as some recent studies 
may imply (Smith et al. 2008; Miller et al. 2009; Quimby et 
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al. 2009), then this may shghtly raise the relative fractions 
of SNe Iln and II-L compared to our study. Recently com- 
missioned untargeted surveys can help constrain this con- 
tribution (see Arcavi et al. 2010, as noted above regarding 
broad-lined SNe Ic in dwarf hosts). 

Our volume-limited survey within 60 Mpc includes 
80 CCSNe, compared to the heterogeneous volume-limited 
study of 92 CCSNe within 28 Mpc summarised by Smartt 
(2009). However, because the LOSS survey was conducted 
with the same telescope in a systematic way, we are able 
to make proper corrections for the observing biases, as Pa- 
per II describes in detail. We also have much more com- 
plete spectroscopic follow-up observations and we monitor 
the photometric evolution of the SNe we discovered, which 
particularly affects the relative fractions of SN II-P vs. II-L, 
Tin, and lib, all of which are sometimes called simply "Type 
II" in initial reports. Thus, samples of SNe using identifica- 
tions from initial reports are often unreliable or unspecific, 
but our study resolves this issue because our more exten- 
sive photometric and spectroscopic follow-up observations 
allow us to more reliably place the SNe in subclasses. Con- 
sequently, our observed fractions of CCSN types differ from 
those of previous studies in a few key respects. The main dif- 
ferences compared to SN fractions listed in various studies 
reviewed by Smartt (2009) are as follows. 

(1) We find a lower SN ILP fraction of only ~48 %, in 
contrast to larger values of 59% in previous studies, although 
some of these did not differentiate among SN II subtypes. 
This impacts the "RSG problem" as discussed below. 

(2) We find correspondingly larger fractions of SNe II- 
L, Iln, and lib compared with Smartt et al. (2009). This 
mostly reflects our spectroscopic and photometric follow-up 
observations mentioned above. 

(3) We find a larger fraction of SNe Ibc than Cappel- 
laro et al. (1997), although similar to other estimates (van 
den Bergh et al. 1987; Prantzos & Boissier 2003; Prieto et 
al. 2008; Boissier & Prantzos 2009; Smartt et al. 2009). The 
number ratio of SNe Ibc to all SNe II that we measure is 
Njhc/Nii = 0.35, whereas Cappellaro et al. found a value 
for Nibc/Nii of only 0.29. Prieto et al. (2008j3 noted that 
Nihc/Nii — 0.27 in the full sample they considered, but they 
also found a metallicity dependence, with higher values com- 
parable to ours at around solar metallicity. The high ratio 
we find is the crux of the "WR problem" that we discuss 
herein. 



3 THE IMF AND PROGENITOR MASSES 

The IMF describes the relative number of stars as a function 
of initial mass, N{m), and within a given mass range this 
dictates the distribution of initial masses for progenitors of 
SNe. We adopt a simple approximation of the IMF as a 
single power law and exponent 7 given by 

N{m) = Cm^, (1) 

where C is a constant. To understand the implications of SN 
rates for massive stars, we investigate the IMF within a mass 

^ Note, however, that Prieto et al. used the Sternberg Astronom- 
ical Institute (SAI) SN catalog, which is a heterogeneous sample 
with unknown systematic biases. 



range bounded by the lowest initial mass that results in a 
CCSN, AfsN, and extending up to the upper mass limit for 
the initial masses of stars. One expects Msn to be around 8 
M0, but there are uncertainties involved, as discussed fur- 
ther below. We take the upper limit to initial masses to be 
150 M0 (Figer 2005), although this choice has little effect 
on our analysis because the most massive stars are so rare 
in the local universe (all the stars from 100 to 150 Mq make 
up less than 2% of the population, comparable to our un- 
certainties). A handy quantity is Fm, which we define as the 
fraction of all CCSNe contributed by stars with initial mass 
m or higher, up to 150 Mq. For an unbroken power-law IMF, 
this is given by 

(■150 

F„ = / N{m')dm (2) 

J m 

[1 - (m/150)^+^l 
[1 - (Msn/150)^+i] ' ^ ' 

where 7 = —2.35 for a standard Salpeter (1955) mass func- 
tion (note that this differs from the logarithmic form that is 
sometimes used, where Salpeter corresponds to F = —1.35). 
Bastian et al. (2010) have provided a recent review of the lit- 
erature on possible variations in the IMF, and conclude that 
there is no clear evidence that the IMF varies strongly in 
the modern universe. Clearly, 1 — Fm is the cumulative frac- 
tion contributed by stars between the lower bound (Msn) 
and m. This assumes that SN progenitors occupy the full 
mass range from Msn to 150 Mq, with no large mass inter- 
val where stars consistently collapse directly to a black hole 
without any visual display (Fryer 1999); the latter remains 
a possibility, and implications are discussed later. 

Figures [5] and [3] show plots of 1 — Fm and Fm , respec- 
tively, for three different representative values of Msn = 8.0, 
8.5, and 9.0 Mq, as well as for two different values of 7 = 
-2.35 (Salpeter 1955) or -2.4 (e.g., Humphreys & McElroy 
1984) for comparison. One can see that small variations in 
7 have little effect on the results. Figure [2] also illustrates 
a hypothetical case of 7 = —1.8, which is large enough to 
make a substantial difference (this is the slope that would be 
needed to reconcile the disagreement between the observed 
fraction of SNe II-P and the observed mass range for the 
corresponding progenitors; see below). This slope, however, 
is more top-heavy than allowed by measurements of local 
stellar populations outside of the inner parts of the densest 
star clusters (see Bastian et al. 2010). 

Small differences in the adopted value of Msn can have 
a substantial effect, however. This is due to the fact that 
lower-mass stars are so much more numerous in a bottom- 
heavy IMF, and small changes in Msn therefore have a 
disproportionate infiuence on the distribution of SN types. 
This is relevant in regard to the still-uncertain lower bound 
to initial masses that experience Fe core collapse and those 
that may suffer less energetic explosions via electron-capture 
SNe (ECSNe). According to Smartt (2009), directly ob- 
served RSG progenitors of normal SNe II-P extend down 
to around 8 Mq and their statistical distribution favours 
MsN = 8.0 ± 1.0 Mq. On the other hand, theories for EC- 
SNe predict that these explosions occur somewhere in the 
range 8-11 Mq depending on assumptions about metallicity, 
mass loss, and other factors (Nomoto 1984; Woosley et al. 
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2002; Kitaura et al. 2006; Wanajo et al. 2009; Pumo et al. 
2009). Theory generally predicts that if ECSNe occur within 
this range, they would tend to be less energetic and fainter 
than a standard Fe CCSN, releasing ^10^" erg of kinetic 
energy (instead of ~10^^ erg) and producing less ^''Ni than 
standard CCSNe. It has been hypothesised that an ECSN 
may have given birth to the Crab Nebula (Davidson et al. 
1982; Nomoto 1984; Kitaura et al. 2006; Wanajo et al. 2009). 

A note of caution is that if the corresponding observed 
visual displays are indeed much fainter than normal CCSNe, 
then it is possible that some of these ECSNe may not be in- 
cluded in the LOSS sample, since we chose to exclude faint 
Type Iln events such as SN impostors, rj Car analogs, LBVs, 
or other peculiar faint transients in the observed fractions of 
CCSNe (Paper II). On the other hand, if ECSNe do not give 
rise to these SN impostors, but appear instead as the rela- 
tively faint end of the distribution of SNe II-P (objects like 
SN 2005cs; Pastorello et al. 2007), then they will be included 
in the LOSS rates as SNe II-P. The luminosity functions in 
Paper II reveal an enhancement at the very bottom of the lu- 
minosity range of SNe II-P. This uncertainty is unfortunate, 
but the ECSN phenomenon is not understood sufficiently 
well to confidently account for it. For this reason, Figures [2] 
and [3] show values for Msn of 8.0, 8.5, and 9.0 Mq and the 
range of uncertainty that this introduces. 

Independent of the questions surrounding ECSN the- 
ory, however, an empirical value of Msn = 8.0±1.0 M0 
is favoured by Smartt (2009) based on the distribution of 
masses for directly observed SN II-P progenitors (although 
one must remember that this value is model dependent as 
well, and subject to systematic effects; see Smartt 2009 for 
details). We adopt Msn = 8.5 Mq for most discussion in this 
work. If the ECSN phenomenon occurs above 8.5 Mq, we 
consider it likely that those ECSNe will be included among 
the population of faint SNe II-P anyway, while those below 
could be excluded if they masquerade as faint transients or 
SN impostors (e.g., Thompson et al. 2009). 

3.1 The RSG Problem 

Red supergiants (RSGs) represent the expected endpoint of 
post-main-sequence stellar evolution for the majority of sin- 
gle stars with initial masses above 8 Mq , and it is straight- 
forward to associate their extended H-rich envelopes with 
SNe II-P — the most common type of CCSN. This has 
long been expected (e.g., Falk & Arnett 1977; Litvinova & 
Nadyozhin 1983; Doggett & Branch 1985; Wheeler & Swartz 
1993), but the RSG/II-P connection has received firm foot- 
ing in the past decade with the identification of RSGs as 
the progenitor stars of several SNe II-P. This work has been 
based on attempting to locate progenitor stars (or upper 
limits to them) in pre-explosion archival data at the same 
position as the SN (Earth et al. 1996; Van Dyk et al. 1999, 
2003a, 2003b, 2003c; Smartt et al. 2001, 2002, 2003, 2004; 
Li et al. 2005, 2006, 2007; Maund & Smartt 2005), and in 
some cases the RSG disappears after the SN has faded. 

These multiple progenitor studies have reassured us 
that RSGs are the progenitors of SNe II-P, but what range 
of initial masses do they imply? There are many potential 
systematic errors involved: masses derived from progenitor 
luminosities rely upon model-dependent evolutionary tracks, 
and circumstellar dust that may have surrounded the pro- 
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Figure 2. The quantity 1 — Fm- This is a cumulative distribution 
function beginning at the bottom of the mass range for CCSNe 
(A/sn), showing the fraction of CCSNe contributed by stars in the 
mass range from MgN up to m. The three solid black curves are 
for three example values of MgN = 8.0, 8.5, and 9.0 Mq using a 
Salpeter slope of 7 = —2.35. The dotted curves are for 7 = —2.4, 
whereas the grey dot-dashed curve illustrates the hypothetical 
top-heavy case of 7 = —1.8 (see text). The long-dashed curve 
labeled "74%" shows the same function for MgN = 8.5 Mq and 
7 = -2.35, but it excludes 26% of the total number (26% is the 
sum of the fractions of all SNe Ibc), assuming that they follow 
a different evolutionary path in close binaries over the full mass 
range considered; this possibility is discussed later in §4.2 and 
4.3. The gray box denotes the range of uncertainty in the upper 
bound to RSG progenitors of SNe II-P, based on the properties of 
progenitors detected so far (Smartt et al. 2008). The horizontal 
line is the observed fraction of SNe II-P. 

genitor could have been vaporised by the SN, causing the ex- 
tinction derived toward the SN progenitor — and therefore 
its luminosity and mass — to be underestimated. Smartt 
(2009) has reviewed the recent literature on the identifica- 
tion of SN II-P progenitors as RSGs in pre-explosion data 
and discussed these systematics. Altogether, Smartt (2009) 
argues that the available collection of SN II-P progenitor 
detections and upper limits favours 8.5-16.5 Mq for the 
range of initial masses, adopting a normal Salpeter IMF, 
and Smartt et al. (2009) give an upper limit to initial masses 
of SN n-P progenitors of 21 Mq with 95% confidence. The 
upper limit in the range of 16.5 to 21 Mq is shown by the 
gray shaded area in Figure [2] 

In our volume-limited sample of SNe, we find that 
SNe n-P constitute about 48% of CCSNe (Figure [g. This is 
a lower fraction than reported in previous studies (Smartt 
2009; Smartt et al. 2009; note that several other previous 
studies did not explicitly separate SNe II-P from II-L or 
other SNe II). Figure [2] compares this LOSS observed frac- 
tion of SNe II-P, Fii-p, to the quantity 1 — Fm (see Eqn. 
3), which is the fraction of CCSNe one expects integrating 
from the bottom of the CCSN range at Msn = 8.0, 8.5, or 
9.0 Mq up to mass m. 

From Figure[5]we see that the initial mass range of 8.5- 
16.5 Mq over which RSG progenitors of SNe II-P have been 
identified would provide more than enough SNe to account 
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for the observed fraction of SNe II- P, under the assumption 
that all stars within this mass range explode as SNe II-P. 
In fact, stars in the initial mass range 8.5-16.5 Mq would 
constitute roughly 62% of all the stars above 8.5 (for 
7 = —2.35) that undergo core collapse, producing too many 
SNe II-P. The mass range 8.5-13.7 would be sufficient 
to produce the observed fraction of SN II-P. 

Thus, there is apparently no RSG problem from the 
"supply-side" point of view, in the sense that the observed 
range of masses for SN II-P progenitors supplies a large 
enough fraction of CCSNe. Looking more closely, however, 
there is a "demand-side" problem in the sense that stars in 
the initial mass range of 8-17 M0 — which are in fact ob- 
served to explode as SNe II-P — produce too many SNe II- 
P compared to the observed fraction of this SN subtype. 
Smartt et al. (2009) did not emphasise this discrepancy in 
their study, presumably because they concluded that SNe II- 
P constitute a larger fraction (~59%) of CCSNe, which 
would be in reasonable agreement with the observed mass 
range within their uncertainties, compared to our value of 
^48%, which is discrepant. A suggestive solution is given 
by the dashed line in Figure [2j which brings the observed 
fraction of SNe II-P and the mass range of detected SN II-P 
progenitors into agreement. This curve is the same as the 
value of 1 — Fm shown by the black curves (7 = —2.35; Msn 
= 8.3 Mq), but multiplied by 74%. The motivation for this is 
that it assumes that the 74% of CCSNe that are Type II are 
distributed evenly across all initial masses, and that there- 
fore the 26% of all SNe that are SNe Ibc have some differ- 
ent origin which is also distributed across all initial masses. 
This might be the case, for example, if all SNe Ibc arise from 
RLOF in binary systems. We return to this question later. 

Note that our comment about the lack of any "supply- 
side" RSG problem is different from the RSC problem 
pointed out by Smartt et al. (2009), which has to do with 
the fact that RSG stars in the initial mass range 17 < m < 
25 Mq are observed to exist, yet they appear to be missing 
from the population of nearby SN II-P progenitors detected 
in pre-explosion data. Smartt et al. (2009) hypothesised that 
these missing progenitor stars may collapse directly to black 
holes without producing successful SNe. Another possible 
solution to this discrepancy, however, is that RSGs in this 
upper mass range continue to evolve into other types of pro- 
genitor stars before core collapse, such as yellow supergiants 
(YSGs), blue supergiants (BSGs), low-luminosity LBVs, or 
Wolf-Rayet (WR) stars, producing SNe of Types II-L, lib, 
Iln, lb, or Ic. Smartt et al. (2009) mentioned this hypothesis 
but disfavoured it, in part because the number of SNe lib + 
II-L -I- Iln was not enough to make up for the missing pop- 
ulation of RSGs, plus other reasons concerning the inferred 
masses for LBVs and progenitors of events like SN 1993J and 
SN 1980K. However, we find that these arguments rely on 
unreliable assumptions and that they provide no compelling 
argument against the idea that RSGs in the initial mass 
range 17-25 M© may continue to evolve before exploding. 

Furthermore, in the volume-limited sample from LOSS, 
we find that SNe II-L, lib, and Iln make up a larger fraction 
of the total SN II group, and SNe Il-P have a lower fraction, 
compared to the study of Smartt et al. (2009). With this 
LOSS sample, we find that there are plenty of remaining 
SNe II besides SNe II-P to account for SNe resulting (even- 
tually) from RSGs known to occupy the higher mass ranges 



above 17 Mq. Another objection stems from the assumption 
by Smartt et al. (2009) that LBVs (the likely progenitors 
of SNe Iln) arise exclusively from stars with initial masses 
above 40 Mq, but there is also a well-known population of 
lower-luminosity LBVs that are thought to be stars in a 
post-RSG phase with initial masses of 20-40 Mq (see Smith 
et al. 2004). In addition to LBVs, Smith et al. (2009) noted 
that the most extreme class of RSGs with high mass-loss 
rates and initial masses of 25-35 Mq could give rise to the 
lower-luminosity SNe Iln. Thus, these considerations alle- 
viate two key objections to the idea that 17-25 Mq stars 
produce other types of SNe that are not Type II-P. 

In fact, there is growing empirical evidence that this is 
indeed the case, supported by direct detections of progenitor 
stars of SNe ILL, lib, and Iln (and ILpec). SN 2009kr is the 
first luminous SN II-L to have a progenitor star identified in 
pre-explosion images (Elias-Rosa et al. 2010b; Fraser et al. 
2010), and it appears to be a YSG. Efias-Rosa et al. (2010b) 
estimate a likely initial mass for the YSG progenitor star of 
18-24 Mq , and infer that it may bridge a gap in progenitor 
mass between SNe II-P and the more massive LBV progen- 
itors of SNe Iln (see below). SN 2009hd in M66 also had a 
Type II-L spectrum, for which Elias-Rosa et al. (2010c) have 
identified another likely YSG progenitor, suggesting an ini- 
tial mass in the range 20-25 Mq. SN 2008cn is yet another 
possible YSG progenitor of a luminous SN II-P (Elias-Rosa 
et al. 2010a), although the large distance to this SN makes 
the progenitor identification less secure. 

Of course, the first SN to have a progenitor identified in 
pre-explosion data was SN 1987A, whose classification was 
Type Il-pec, and which was inferred to have an ~18 Mq 
BSC progenitor that was in a post-RSG phase (see Arnett 
1987; Arnett et al. 1989). The progenitor of the SN lib 1993J 
was inferred to be a Mq « 15 Mq K-type RSG with a 
large radius but small H envelope mass (Aldering et al. 1994; 
Filippenko et al. 1994; Van Dyk et al. 2002; Maund et al. 
2004). In both cases, binary evolution was invoked to explain 
the status of the progenitors at the time of core collapse 
(Nomoto et al. 1993; Podsiadlowski et al. 1993; Aldering et 
al. 1994; Woosley et al. 1994). Lastly, so far only one SN lln 
(SN 2005gl) has a progenitor star identified in pre-explosion 
data, and it was a massive LBV (Gal- Yam & Leonard 2009). 

Collectively, these results argue that RSGs in the range 
of masses above the observed range for SN Il-P progeni- 
tors may indeed continue to evolve after the RSG phase due 
to further mass loss (in either single- or binary-star evolu- 
tion) , to produce other types of SNe. This relieves the RSG 
problem proposed by Smartt et al. (2009), and removes the 
empirical motivation for inferring that massive stars in some 
mass range collapse directly to a black hole (BH) without a 
visible SN display. In fact, we find that the latter inference 
would introduce other problems that are at odds with the 
observed fractions of CCSNe, as discussed further below. 



3.2 The WR Problem 

Unlike the detected progenitors of SNe II-P, II-L, lib, and 
lln, the progenitor stars that have shed their H envelopes 
to make SNe Ibc have never yet been identified directly in 
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Figure 3. The quantity Fm (Eqn. 3). This is similar to Figure[2l 
but with a cumulative distribution funetion beginning at the up- 
per mass limit and working down, showing the fraction of stars 
in the mass range between m and the assumed upper mass limit 
at 150 Mq. The horizontal lines mark the observed fractions of 
SNe Ic and the sum of SNe lb + Ic + Ibc-pec. The vertical dashed 
line at Mwr = 35 marks the initial mass above which H-free 
WR stars are thought to originate, inferred from observations. 
The grey curve is the same as the solid black curve corresponding 
to MgN = 8.5 Mq, but multiplied by 0.26 to mimic the distri- 
bution of SNe Ibc if they were evenly distributed across the full 
mass range. 

pre-explosion data0 The known stars that most naturally fit 
the bill for progenitors of SNe lb and Ic are the Wolf-Rayet 
(WR) stars of the WN and WC subclasses, respectively, be- 
cause of their relatively H-free surface composition and their 
small stellar radii (e.g., Woosley & Bloom 2006). The distri- 
bution of WR stars in galaxies appears marginally consistent 
with that of SNe Ibc, although this depends on metallicity 
(Leloudas et al. 2010)0 

It is straightforward to expect that the WC subclass 
would explode to produce SNe Ic, but it is not so clear if the 
WN subclass explodes as SNe lb, or if instead the WN stars 
should continue to evolve by virtue of their own mass loss to 
become WC stars before exploding as SNe Ic. Evidence that 
some WN evolve to WC is that the WN/WC ratio is ~1 
in the Milky Way and higher at lower metallicity, whereas 
from LOSS we find that Nib/Nic is only < 0.5. This interplay 
may be luminosity and metallicity dependent (as discussed 
further below), and comparisons of WR and SN Ibc positions 
in galaxies give mixed results (Leloudas et al. 2010). The fact 

^ One putative exception is SN 2010O, for which Nelemans et al. 
(2010) detect a variable X-ray source at the SN position in pre- 
explosion data. Nelemans et al. (2010) claim that this may have 
been a WR/black hole binary system, where SN 2010O was the 
second SN in the system when a WN star produced the observed 
SN lb. The WN star itself was not detected, however. 
^ We do not include the group of luminous H-rich late-type WN 
stars, or WNH stars (Smith & Conti 2008), which are probably 
still in core-H burning and are more like O-type stars with en- 
hanced winds. 



that no normal SN lb or Ic has an identified progenitor staiQ 
makes the identification of luminous WR stars as the only 
progenitors of both SNe lb and Ic uncertain. 

Standard single-star evolution models (e.g., Meynet et 
al. 1994; Heger et al. 2003) predict that strong line-driven 
stellar winds at high luminosity will cause stars more mas- 
sive than some threshold mass to completely shed their H 
envelopes. This leaves He cores that are observed as WR 
stars (e.g., Conti 1976), and which should explode to make 
SNe Ibc. For convenience, we define Mibc as the initial mass 
dictated by the observed fraction of SNe Ibc, above which all 
progenitors have fully shed their H envelopes before core col- 
lapse. From Figure O we find that SN statistics from LOSS 
show that Afibc ~ 22M0 . 

Similarly, we define Mwa as the initial stellar mass 
above which a massive star is expected to shed its H en- 
velope. If standard single-star evolution applies, then we 
should find Mibc = Mwr- However, standard single-star 
evolutionary models such as those by Heger et al. (2003) 
predict a much higher value of Mwr = 34 Mq at solar metal- 
licity, and they suggest that Mwr rises to even higher initial 
masses at lower metallicity due to the strong metallicity de- 
pendence of line-driven winds that are assumed to dominate. 
A problem recognised in recent years is that these single- 
star evolutionary models have used empirical prescriptions 
for mass-loss rates that are now known to be far too high 
by factors of 3-10 compared to observed mass-loss rates, 
as noted in §1. Using more realistic wind mass-loss rates 
would change the predictions significantly, such that single- 
star evolution would not be able to account for the popu- 
lation of WR stars or SNe Ibc, even at solar metallicitylfl 
Smith & Owocki (2006) have discussed this, pointing out 
that giant LBV-like eruptions may provide a way to make 
up the deficit, but the mass range over which this applies is 
uncertain; eruptive mass loss is probably dominant in only 
the most massive stars. Observations of WR stars associated 
with star clusters suggest a value of AIwb. around roughly 
35-40 Mq for most WR stars in the Milky Way (Schilde 
& Meader 1984; Massey et al. 1995, 2001; Crowther et al. 
2006; Crowther 2007; Massey 2003; Humphreys et al. 1985). 
We therefore adopt Mwr ~ 35 Mq for the majority of WR 
stars at roughly solar metallicity|j 

A pre-explosion source was identified at the position of the fa- 
mous object SN 2006jc, but this is a highly unusual case. The 
pre-explosion object was seen only as a transient source in a brief 
eruptive phase 2 yr before the SN, and the subsequent SN was 
a very unusual event, probably an underlying SN Ic whose shock 
overtook a dense He-rich circumstellar shell to produce a SN Ibn 
with bright, relatively narrow He I lines in the spectrum (see 
Pastorello et al. 2007; Foley et al. 2007; Smith et al. 2008a). Nev- 
ertheless, the ~10^ km s"'^ speed of the pre-shock CSM seems 
consistent with a compact WR-like progenitor. 
^ Models with rotation (e.g., Meynet et al. 2008) have been pro- 
posed to yield lower values of Mwr as low as ~25 Mq that are in 
better agreement with Mjbc, but this is because they have even 
higher mass-loss rates, violating observational constraints even 
more severely; this makes rotation an unlikely solution. 
^ There is also evidence that some lower-luminosity early-type 
WN stars may originate from lower initial masses down to ~25 
Mq (Crowther 2007), but stellar winds at solar metallicity (and 
probably even LBV eruptions as well) are insufficient to strip 
their H envelopes. In close binary evolution, however, complete 
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Figure |3] highlights a serious problem with assigning 
classical WR stars as the exclusive progenitors of SNe Ibc. 
Namely, the fraction of all stars experiencing CCSN above 
MwR = 35 M© only accounts for about half the number 
needed for the observed fraction of SNe Ibc. To account for 
all SNe Ibc in this simple prescription — where more lumi- 
nous stars have higher mass-loss rates and therefore become 
WR stars and SNe Ibc by virtue of their own mass loss — 
would require that SN Ibc progenitor stars extend from the 
upper mass limit down to around 22 Mq. In other words, 
the WR problem can be stated simply as 

MwR >> Mibc. (4) 

According to Figure [S] roughly half the SN Ibc popu- 
lation must originate from stars that are less massive than 
initial masses corresponding to the observed or theoretically 
expected population of WR stars. There are several possible 
solutions to this problem. (1) The WR phase for many lower- 
mass stars is not observed because it is extremely short- 
lived, perhaps because eruptions in late evolutionary phases 
remove the remaining H envelope even down to lower masses 
than we normally associate with LBVs. In this case, how- 
ever, some SNe Ibc should show signs of interaction with 
H-rich CSM at late times, because that H must have been 
shed recently, while only a few do. SN 2001em is one ex- 
ample (Chugai & Chevalier 2006; Schinzel et al. 2009; Van 
Dyk et al. 2009), but perhaps there are more where the 
CSM interaction is missed at very late times. (2) Alterna- 
tively, the population of H-free stars that correspond to the 
progenitors of almost half of SNe Ibc may be underlumi- 
nous because of significant mass loss in binary RLOF. If 
under luminous, their radiation-driven winds — and hence, 
their emission-line spectra — may be weak and so they are 
not discovered or identified observationally as classical WR 
stars. These may be hidden by brighter companion stars in 
binary systems (i.e. the overluminous mass gainers), mak- 
ing them more difficult to observe. Smartt (2009) mentioned 
this as a potential explanation for the lack of any detection 
of SN Ibc progenitors so far. The idea that RLOF dominates 
the population of SN Ibc progenitors was suggested long ago 
(e.g., Filippenko 1991; Branch, Nomoto, & Filippenko 1991; 
Podsiadlowski et al. 1992) but has been hard to confirm. We 
find, as discussed further below, that this is the likely ongin 
of at least half and possibly most SNe Ibc. 

This is not to say that the expectations of single-star 
evolution are completely irrelevant. While binary RLOF 
may be largely independent of metallicity and initial mass, 
stellar winds may still play an important role. RLOF in bi- 
naries provides a likely way to strip the H envelopes at any 
metallicity, but it is less likely to strip the He envelope except 
for the shortest-period systems. The same goes for shedding 
the H envelope via giant LBV-like eruptions, which may 
also be insensitive to metallicity (Smith & Owocki 2006). 
However, the subsequent evolution of the stripped He core 
— from one with a small residual H mass to a H-free and 
He-rich surface, and eventually toward removal of the He 
layer as well — can be accomplished by the line-driven wind 



removal of the H envelope can occur at much lower initial masses 
down to ~15 Mq or less (e.g., Eldridge et al. 2008; Podsiadlowski 
et al. 1992). We will return to this later. 



of the WR star itself, which does depend strongly on both 
metallicity and luminosity (and therefore initial mass). This 
is supported by the observation that the ratio of WN stars 
to WC stars varies from ~1 in the Milky Way to 5 and 10 
in the LMC and SMC, respectively (e.g., Crowther 2009). 
Even as binaries, some of the SNe Ibc - in particular the 
SNe Ic - may therefore appear to obey expected trends of 
single-star evolution, where the most luminous and higher 
metallicity stars are more able to shed their He envelopes via 
radiation-driven winds or eruptions, leading to WC stars and 
SNe Ic. This may explain why studies of the positions of WR 
stars and SNe in their host galaxies find that WC stars and 
SNe Ic seem to imply higher initial mass and higher metal- 
licity environments (Kelly et al. 2008; Anderson & James 
2009; Leloudas et al. 2010; see also Papers I, II, and HI), 
even if binary evolution or LBV eruptions dominate the re- 
moval of the H envelope. We emphasise that it will be quite 
important in future studies to distinguish between SNe lb 
and Ic while studying SN statistics as functions of metallic- 
ity and redshift. By the same token, it may be important to 
distinguish among subtypes of SNe II, as discussed next. 

3.3 The LBV Problem 

At odds with the standard scenario for the formation of 
WR stars as the descendants of the most massive stars with 
Mo ^ 35 Mq is the uncertain fate of LBVs and their connec- 
tions to SNe. In this standard scenario (e.g., Conti 1976), 
winds of O-type stars on the main sequence shed much of 
the H envelope, leaving a very brief ('^lO* yr) transitional 
LBV stage at the end of core-H burning that finishes the job 
of forming H-free WR stars. Recent work (Bouret et al. 2005; 
FuUerton et al. 2006; Puis et al. 2006) has demonstrated that 
O-star winds are clumped and that their mass-loss rates are 
too weak, so it appears likely that LBV giant eruptions must 
dominate this mass loss if WR stars are to form via single- 
star evolution (Smith & Owocki 2006). If these eruptions are 
not strong enough, the star will fail to shed much of its H 
envelope before core collapse, producing a SN Iln (Smith & 
Owocki 2006). 

In fact, recent studies have provided mounting evidence 
that some LBVs explode as SNe before the stars are able 
to fully shed their H envelopes. Luminous SNe Iln, which 
are thought to be powered by shock interaction with dense 
CSM, require large masses of material ejected in sudden 
eruptions that occur within decades before core collapse, in 
some cases as high as 10-20 Mq (Chugai et al. 2004; Smith et 
al. 2007, 2008b, 2010; Woosley et al. 2007; Smith 2008). The 
large CSM masses for luminous events like SN 2006tf and 
SN 2006gy require very massive progenitor stars to account 
for the mass budget, since the large ejecta mass corresponds 
only to the H-rich envelope ejected just before core collapse 
(i.e., the true initial mass of the star also includes the He 
core and any mass shed during the star's lifetime). One of 
the hypotheses for the pre-SN mass ejections of SNe Iln 
is that they suffered pulsational pair instability ejections 
before core collapse, in which case very massive progenitor 
stars with Mq > 95 Mq are needed (Heger et al. 2003) . 

There are other, anecdotal signs of a link between LBVs 
and SNe Iln as well, having to do with their wind speeds, ab- 
sorption profiles, and circumstellar nebulae (Kotak & Vink 
2006; Smith 2007; Trundle et al. 2008). Much more directly. 
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Figure 4. Mass ranges implied by the observed fractions of SN types in a standard single-star evolutionary framework, where higher 
initial masses lead to higher mass-loss rates, and consequently greater stripping of the H and He envelopes. Afibc ~ 22 is the dividing 
point above which stars must fully shed their H envelopes in this scenario. 



Gal- Yam & Leonard (2009) showed that the LBV-like pro- 
genitor of the SN Iln 2005gl subsequently disappeared, pro- 
viding a strong case that LBVs do in some cases explode 
as SNe Iln, despite that fact that no contemporary stel- 
lar evolution models predict this. Gal- Yam & Leonard in- 
ferred a high initial mass of > 50 Mq for the progenitor of 
SN 2005gl. The "LBV problem," then, is the fact that LBVs 
or some other very massive, unstable H-rich stars explode 
as SNe Iln, even though current models expect very massive 
stars to shed their H envelopes. 

If SNe Iln truly arise from massive LBV progenitors, 
exactly what ranges of initial mass are required? How can 
we divide the IMF such that very massive progenitor stars 
can yield both SNe Iln and the SNe Ic that are supposed 
to come from the WC descendants of very massive stars? 
What scenarios are consistent with the observed fractions of 
various types of CCSNe? We investigate this problem next. 



4 HYPOTHETICAL SCENARIOS FOR 
DIVIDING THE IMF 

Given the problems and complications between progenitor 
scenarios expressed in the previous section, we now address 
the problem from a simpler empirical point of view. Here 
we ask how one can subdivide the IMF of massive stars in 
a way that is consistent with the fractions of various CCSN 
types observed in LOSS, while also meeting requirements 
imposed by our knowledge of the likely progenitor stars. For 
simplicity, in all cases we adopt a Salpeter IMF within the 
mass range bounded by the lowest initial mass for which SNe 
occur, assumed to be Msn = 8.5 Mq, up to the proposed 
upper mass limit for initial masses at 150 Mq (Figer 2005). 

This is meant to be exploratory and demonstrative, 
rather than definitive. We consider extreme hypotheses such 



as one where all massive stars obey expectations of single- 
star evolution (e.g., Heger et al. 2003), and alternatively, 
where all stripped-envelope SNe arise from binary RLOF 
(e.g., Filippenko 1991; Podsiadlowski et al. 1992), and we 
evaluate merits and drawbacks of each. We also mention a 
compromise "hybrid" scenario. Our analysis is intended to 
guide intuition in future studies, and to provide tests for 
single/binary-star population synthesis models. 



4.1 Dominated by Standard Single-Star Evolution 

We first consider the familiar hypothesis that at a given 
metallicity, increasingly more massive and more luminous 
stars have monotonically increasing mass-loss rates, such 
that higher initial masses invariably lead to greater strip- 
ping of the H and He envelopes. It is essentially a hypoth- 
esis that single-star mass loss dominates over close binary 
interactions in stripping a massive star's envelope, thereby 
determining the distribution of SN types. This is widely con- 
sidered to be the "standard" view of mass loss connecting 
stellar initial masses to their ultimate fates as a function 
of metallicity (e.g., Heger et al. 2003). In this picture, the 
most massive stars fully shed their H envelopes by virtue 
of their own strong winds or LBV-like eruptions to produce 
SNe Ibc. At intermediate masses, stars do not fully shed 
their H envelopes, instead producing SNe Iln, lib, and II-L, 
depending on how much H mass was lost, and how recently 
this occurred (i.e., the density of the immediate CSM). The 
lowest mass range corresponds to RSGs that do not shed 
their H envelopes and produce SNe II-P. 

Figure |4] shows how the IMF would need to be subdi- 
vided in this hypothetical single-star framework, dictated by 
the observed fractions of various SNe types determined by 
LOSS (Figure[T]). Figure|4]is largely a more succinct restate- 
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ment, in graphical form, of the discussion above concerning 
inconsistencies with RSGs, WR stars, and LBV progenitors. 

While in principle this scenario is consistent with 
the qualitative expectation that more massive stars have 
stronger pre-SN mass loss, it also comes with many inconsis- 
tencies, and conflicts with several observational constraints 
on the likely progenitors of various SN types. Some obvious 
problems evident from Figure [4] are the following: 

(1) The mass range occupied by the observed fraction of 
SNe II-P (8.5-13.7 M©) is too small compared to the directly 
observed mass range of SNe II-P progenitors, 8.5-17 Mq or 
more (Smartt 2009). In other words, if all stars in the range 
8.5-17 Mq produced SNe II-P, then the fraction of CCSNe 
that are II-P would be much higher than observed. 

(2) This scenario contradicts the observational indica- 
tion that some SNe Iln have very massive LBV-like progen- 
itors, as discussed above. A few SNe Iln appear to have pro- 
genitors with initial masses of at least 50-80 M©, whereas 
Figure |4] requires that no stars above 22 Mq retain any H 
envelopes at core collapse. 

(3) Similarly, recent identifications of yellow supergiants 
as SNe II-L progenitors place them at the upper extreme of 
the range allowed for SNe II, or even above 22 Mq. Masses 
inferred for the SNe II-L 2009kr and 2009hd are 18-24 Mq 
and 20-25 Mq, respectively (Ehas-Rosa et al. 2010b; 2010c). 

(4) Most importantly, there are far too many observed 
SNe Ibc, requiring that all stars above 22 Mq completely 
shed their H envelopes, whereas the expected value Mwr is 
roughly 35 Mq. Even at solar metallicity, stars below Mwr 
do not have sufficiently high mass-loss rates to shed their 
H envelopes — certainly not through metallicity line-driven 
winds or RSG winds, and probably not through LBV erup- 
tions either. The known initial mass range for most nearby 
WR stars accounts for only half the SN Ibc population. 

We regard this disagreement as strong evidence that 
standard single-star evolution with mass loss simply cannot 
account for the observed distribution of SNe types, and that 
binary RLOF is therefore needed to account for at least half 
of the SN Ib/Ic population, possibly most of it (this is the 
next case discussed below). 

Introducing the hypothesis that some stars collapse di- 
rectly to BHs without making a visible SN does not help. If 
we assume that the most massive stars collapse to BHs (with 
>40 Mq, for example; Fryer 1999), then it pushes the divid- 
ing mass between SNe II and SNe Ibc to even lower values, 
making the problem worse. It would also push the upper 
mass for SN II-P progenitors even lower, causing an even 
worse discrepancy with direct progenitor mass estimates. A 
partial solution relying on BHs would require a finely tuned 
or carefully chosen set of intermediate mass ranges for BHs, 
but it is still unsatisfying (i.e., assuming that stars of, say, 
20-30 Mq initial mass collapse to BHs could bring the mass 
range of SNe Ibc into better agreement with Mwr, but it 
would worsen the problem in points 1-3 above). Direct SN- 
less collapse to a BH may nevertheless be a possibility. Bet- 
ter constraints on the disappearance of stars without SNe 
are needed (e.g., Kochanek et al. 2008). 

4.2 Dominated by Close Binaries #1 

An alternative to single-star mass loss is that mass ejection 
or mass transfer via RLOF in interacting binaries plays a 



dominant role in stripping away the H envelope for a signif- 
icant fraction of SN progenitors. This binary hypothesis for 
explaining WR stars and SNe Ibc has been around longer 
(Paczyhski 1967) than the idea that stellar winds of single 
stars remove the H envelope (Conti 1976). Several studies of 
the effects of binary RLOF on massive star evolution have 
been conducted (e.g., Podsiadlowski et al. 1992; Wellstein & 
Langer 1999; Vanbeveren et al. 2007; Eldridge et al. 2008). 
It has been difficult to confirm or refute the idea that bi- 
nary RLOF dominates the removal of the H envelopes in 
massive stars because of uncertainties in the binary fraction 
as a function of initial mass (see Kobulnicky & Fryer 2007) 
and the large number of free parameters in binary models. 
Also, until very recently (when mass-loss rates of hot stars 
have been revised downward), single-star evolution seemed 
to provide a sufficiently plausible alternative. We argue here 
that low mass-loss rates of single stars combined with the 
large SN Ibc fraction now demand that binary RLOF plays 
a dominant role for a large fraction of SNe Ibc. 

Figure[5]shows a simplified scenario that is radically dif- 
ferent from Figure 2] It represents the other extreme where, 
instead of assuming that all stars shed their H envelopes via 
their own winds in single-star evolution, we adopt the op- 
posite premise that all SNe Ibc have lost their H envelopes 
via RLOF in binary systems (following Kobulnicky & Fryer 
2007; Fryer et al. 1998, 1999; Filippenko 1991; Podsiadlowski 
et al. 1992; Eldridge et al. 2008). To create Figure[Sl we sim- 
ply assumed that the observed fraction of SNe Ibc, ^26%, 
is identical to the fraction of massive stars that lose their H 
envelopes in RLOF, and that the remaining H-bearing SNe 
are distributed across the full mass range. We of course do 
not know the binary frequency as a function of initial mass, 
so for simplicity, this 26% is then distributed evenly across 
all initial masses of SN progenitors. 

This simple "Binary #1" scenario has some advantages 
over the standard single-star hypothesis, as well as some 
drawbacks, a follows. 

(1) The initial mass range of 8.5-18 Mq occupied by the 
~48% of CCSNe that are SNe II-P is now in much better 
agreement with the inferred mass range of RSG progenitors 
(Smartt 2009). 

(2) The mass range of SNe II-L is in better agreement 
with recent detections of progenitors mentioned above, al- 
though perhaps somewhat too high, and it is unclear how 
the difference between SNe II-L and lib arises naturally in 
this scenario. 

(3) By redistributing the remaining SNe II over all ini- 
tial masses, this scenario allows for SNe Iln to be associated 
with the most massive stars, consistent with their presumed 
massive LBV progenitors, and with the pulsational pair in- 
stability in the most extreme cases. This scenario also has 
the appealing characteristic that the H-rich sequence II-P 
— >■ II-L/IIb — > Iln corresponds to a sequence of single pro- 
genitors with increasing mass-loss rate, and hence SNe with 
increasing CSM interaction. 

(4) Owing to the fact that SNe lb are relatively rare, the 
mass range of SNe lb (including SNe Ibc-pec) is surprisingly 
narrow and low, at only 8.5-12.4 Mq, if they occupy lower 
masses than SNe Ic within the binary zone in Figure[5l There 
are currently no direct detections of SN lb progenitors. 

A potential objection to this simple binary scenario is 
that the fraction of SNe Ibc compared to SNe II is observed 
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Figure 5. Same as Figure |4] but now assuming that close binary evolution and RLOF is a necessary ingredient to explain the loss of 
the H envelope for all SNe Ibc. In this "Binary #1" scenario, the fraction of all massive stars that lose their H envelopes in this way is 
dictated to be 26%, which is the same as the observed fraction of SNe Ibc. For simplicity, these binaries are divided equally among all 
initial masses; consequently, the remaining stars that fail to shed their H envelopes (all SNe II) are redistributed across the full range of 
initial masses as well, following expectations that more massive stars have higher mass-loss rates. 




Figure 6. Same as Figure O except that we have included SNe lib in the same group with SNe Ibc, all of which are assumed to have 
their envelope stripping dominated by RLOF in close binary systems. 



to be metallicity dependent (Prantzos & Boissier 2003; Pri- 
eto et al. 2008; Boissier & Prantzos 2009; Papers I, II, and 
III), as are different WR subtypes (Crowther 2007), and 
tiiat SNe Ic are thought to be associated with massive stars 
and higher metalhcity because of their specific locations in 
galaxies (Kelly et al. 2008; Anderson & James 2009; Papers 
I and II). However, this may still be true even in the simple 



binary scenario. Binary RLOF is only effective at remov- 
ing the H envelope in most cases, leaving the He core ex- 
posed. More massive and luminous stars will have stronger 
winds with higher radiation-driven mass-loss rates, which 
will dominate the subsequent evolution. Consequently, only 
more massive stars (or perhaps the closest binaries) expe- 
rience further significant mass stripping, driving evolution 
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from WN to WC to produce SNe Ic rather than SNe lb. A 
corollary is that line-driven winds of WR stars are metallic- 
ity dependent, so while removal of the H envelope (either by 
RLOF or LBV-type eruptions) is insensitive to metallicity, 
the further evolution from WN to WC (and hence, the pro- 
duction of SNe Ic) will be highly dependent on metallicity. 

One last complication is that this scenario places 
SNe lib as single stars and SNe lb as binaries. This leaves 
us without a satisfying explanation as to why such a tiny 
difference in surface H mass separates SNe lib and SNe lb 
(e.g., Chornock et al. 2010; Elmhamdi et al. 2006; Filippenko 
et al. 1994), which otherwise look extremely similar, and it 
ignores observational results suggesting that the progenitor 
of the nearby Type lib SN 1993J was most likely a binary 
system (Aldering et al. 1994; Maund et al. 2004; Maund & 
Smartt 2009). There is also evidence for binarity in the case 
of SN 2001ig (Ryder et al. 2004, 2006; Silverman et al. 2009; 
Maund et al. 2007), which was also Type lib. These issues 
motivate the alternative binary scenario discussed next. 

4.3 Dominated by Close Binaries #2 

The second binary-dominated progenitor scenario that we 
consider is similar to the first, except that we now include 
all SNe lib along with SNe lb and Ic as stars that lose their H 
envelopes primarily through binary RLOF. The motivation 
for this, as explained above, is the close morphological rela- 
tionship between SNe lib and lb — SNe lib essentially are 
Type lb except for a small amount of H at early times — plus 
the observational evidence of the progenitor of SN 1993 J and 
models for its evolution that are suggestive of a binary sys- 
tem (Podsiadlowski et al. 1993; Aldering et al. 1994; Maund 
et al. 2004; Maund & Smartt 2009). 

Figure |S] shows how the IMF could be divided accord- 
ing to observed SN fractions if we assume that all "stripped- 
envelope SNe," now including SNe lib along with SNe Ibc, 
arise from binary RLOF. Including SNe lib as binary sys- 
tems has three main consequences compared to the Binary 
#1 scenario: 

(1) The fraction of all CCSNe progenitors that lose their 
H envelopes through binary RLOF is higher, at ~37% in- 
stead of ~26%. Note that both cases Binary #1 and #2 
imply rather high binary fractions, as the stripped-envelope 
progenitors are mainly the mass losers in RLOF binary sys- 
tems, but the implied close binary fraction is within reason 
(see Kobulnicky & Fryer 2008). 

(2) The upper mass bound for SN II-P progenitors is 
shifted to higher masses (23.6 Mq). This upper bound is 
somewhat troublesome, as it exceeds the 95% confidence 
upper limit of 21 M© derived from the properties of SN II-P 
progenitors (Smartt et al. 2009). 

(3) Most significantly, the mass range for SNe lb shifts 
to higher mass progenitors than in the Binary #1 scenario. 
Assuming that progenitors of SNe lb are less massive than 
SNe Ic in the Binary #1 scenario would dictate that SNe lb 
arise from initial masses of 8.5-12.4 M©; as noted above, 
this is low and quite narrow. If we assume the same for the 
Binary #2 hypothesis, but also add the assumption that 
SNe lib, in turn, are less massive than SNe lb, then the 
corresponding ranges of initial masses would be 8.5-11 M© 
for SNe lib, 11-16 Mq for SNe lb, and >16 M© for SNe Ic. 
This is an improvement over the Binary #1 scenario in that 



it pushes the dividing mass between SNe lb and Ic to higher 
masses, although 16 M© still seems quite low for WR stars 
that we expect to shed their own He envelopes via line-driven 
winds. This is remedied in the "hybrid" scenario discussed 
next. An important caveat is that the monotonic transition 
SNe lib — !■ lb — >■ Ic with increasing initial mass is probably 
not strict, as it also depends on initial binary separation 
(i.e., very close binaries can remove all of the H and even 
He layers in RLOF). Thus, SNe lib could extend to higher 
masses than 11 M© if they arise in relatively wide binaries, 
for example (see below). 

By dividing SN types into two different and distinct 
channels corresponding to single stars and binaries, the Bi- 
nary #2 hypothesis has the appealing quality that it pro- 
vides a natural continuity in SN types within each channel, 
which is lacking otherwise. With increasing levels of enve- 
lope stripping due to RLOF followed by WR wind mass loss, 
the binary channel gives SNe lib — ^ lb ^ Ic. There may be 
a continuum of SN progenitors with different levels of enve- 
lope stripping, probably corresponding to increasing initial 
metallicity or luminosity, such that a small amount of resid- 
ual H separates SNe lib and lb (e.g., Elmhamdi et al. 2006; 
Chornock et al. 2010), whereas a small difference in He mass 
may separate SNe lb from Ic. 

In the Binary #2 hypothesis, there is also now a natural 
continuity in the single-star channel, giving SNe II-P — !> II- 
L — !■ Iln with increasing initial mass and pre-SN mass loss, 
and without the puzzling ambiguity between the origins of 
SNe lib and II-L. The few direct detections of progenitors 
that are available support the notion that the progenitors 
of SNe II-L are more massive than those of SNe II-P (Elias- 
Rosa et al. 2010b; 2010c), and that progenitors of SNe Iln 
are more massive than SNe II-L (Gal- Yam & Leonard 2009). 
The same is true for levels of CSM interaction: SNe II-P 
tend to have extremely weak or undetectable CSM interac- 
tion signatures, SNe II-L tend to have stronger radio and 
X-ray emission (Sramek & Weiler 1990), and their Ha pro- 
files with weak P-Cygni features are thought to arise from 
heating of the SN ejecta by CSM interaction (e.g., Chugai 
1991). SNe Iln obviously have the strongest levels of CSM 
interaction, but there is wide diversity even among the sub- 
class, with the faintest SNe Iln like SN 2005ip looking basi- 
cally like a SN II-L with strong narrow emission lines (Smith 
et al. 2009b), whereas the CSM is opaque and qualitatively 
changes the SN in more luminous SNe Iln such as SN 2006tf 
and SN 2006gy (Smith et al. 2008b). The fuU range for 
SNe Iln (34-150 Mq) encompasses the most luminous RSG 
that may be responsible for the fainter SNe Iln (Smith et al. 
2009a; see also Yoon & Cantiello 2010), intermediate cases of 
SNe Iln consistent with normal LBVs (Gal- Yam & Leonard 
2009), as well as the most massive stars with violent pre-SN 
mass loss (Smith et al. 2007, 2010; Woosley et al. 2007). 

If we relax the requirement that all of the most massive 
single stars make successful SNe Iln, then the lower-right 
corner of Figure [6] provides an attractive parameter space 
for massive stars that can collapse to a BH without making 
a SN display. If, for example, we allow all single stars above 
50 Mq in the Binary #2 scenario to quietly make BHs, then 
the redistribution of the remaining mass ranges for SNe II-P, 
II-L, and Iln are still in rough agreement with observational 
constraints. Of course, this would fail to produce the very 
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Figure 7. Our favoured scenario, combining single and binary star evolution. This is the same as Figure |6] except that now we have 
taken roughly half of the SNe Ic (8.8% of all CCSNe, to match the fraction of SNe Iln) away from the binary RLOF population and 
mixed them with the single-star population. SNe Ic that arise from single stars are below the dashed line. Thus, in this scenario we 
assume that half of all single stars above ~23 are able to shed their H envelopes via winds or LBV eruptions, while the other half 
retain their H envelopes until just before core collapse, producing SNe Iln. The difference among the most massive stars would depend on 
the efficiency of winds and LBV eruptions, which in turn may depend on properties such as metallicity or rotation. The specific numbers 
shown here are meant to provide just one example of a potential hybrid scenario. 



luminous SNe Iln that are thought to come from the most 
massive stars. 

A drawback of this Binary #2 scenario is that the initial 
mass range for SNe Ic still reaches to uncomfortably low 
masses, and therefore dominates most of the mass range 
for binary progenitors. Note that if we allow some of the 
most massive stars in the binary channel to undergo quiet 
BH collapse, we would need to shift the boundary between 
SNe Ic and lb to even lower initial masses, exacerbating this 
problem. Also, the Binary #2 scenario does not allow any 
SNe Ic to come from single stars, and begs the question of 
the origin and fate of single WR stars, which presumably 
arise from eruptive LBV mass loss in very massive stars or 
perhaps through strong winds at super-solar metallicity. The 
next scenario allows some of the most massive single stars 
to produce SNe Ic as well. 

4.4 A Hybrid Scenario 

One can, of course, play this game ad nauseum by adjust- 
ing the fraction of SN progenitors that experience binary 
RLOF, and redistributing the remainder among single stars 
in various ways. Figure [7] shows an example of one "hybrid" 
scenario, which is a compromise between the standard view 
of single-star evolution and the Binary #2 scenario. Here we 
have assumed that roughly half of the SN Ic population (we 
take a fraction equal to 8.8% of all CCSNe for convenience, 
equal to the SNe Iln fraction) may arise from single-star evo- 
lution, while the remainder of SNe Ic form via binary RLOF 
along with SNe lb (including Ibc-pec) and SNe lib as before, 
so that the binary RLOF fraction is 28% in this hypotheti- 



cal scenario. The binary fraction may be somewhat different 
or may be mass dependent, and one can adjust a version of 
Figure [7] accordingly to match precise values; the goal here 
is to be conceptual. 

Although such a scenario may seem more complicated 
and somewhat ad hoc, it is well motivated, and balances 
several competing factors. Among the most massive stars 
with initial masses above 23 M0, it allows single stars to 
die as either SNe Ic or Iln. This may be the case if the effi- 
ciency of single-star mass loss depends on additional factors 
such as rotation or metallicity. One can imagine, for exam- 
ple, that very massive stars may be unable to shed their 
H envelopes if low metallicity or slower initial rotation rates 
weaken their winds or tame the LBV instability. Under these 
circumstances, massive stars might then die as SNe Iln if 
they suffer core collapse while still in the process of attempt- 
ing to shed their H envelopes. Indeed, we noted in Paper II 
that SNe Iln tend to prefer smaller, lower metallicity galax- 
ies. The remainder of more rapidly rotating single stars or 
higher metallicity single stars might successfully shed their 
H envelopes via winds or LBV eruptions and die as SNe Ic. 
LBV eruptions do seem to be more catastrophic among the 
most massive stars (Smith & Owocki 2006). 

Aside from being hypothetical, this scenario has no ob- 
vious disadvantages in view of our knowledge of SN progen- 
itors, and it has some strengths as follows. 

(1) It maintains very good agreement between the mass 
range of SNe II-P and the inferred mass range of directly 
detected RSG progenitors (Smartt 2009). Putting some of 
the SNe Ic back into the single-star channel has the con- 
sequence that it lowers the upper mass bound required for 
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Figure 8. This is virtually the same as Figure[7] except that the SNe Iln and single-star SNe Ic are not divided equally across the range 
of masses; instead, the SNe Iln occupy lower masses than single-star SNe Ic. For equal fractions, the dividing mass between SNe Iln and 
single-star SNe Ic would need to be ~ 37 Mq. 



SNe II-P compared to the Binary #2 scenario, improving 
the agreement with observations. Obviously, we could have 
chosen the fraction of SNe Ic to be a little larger in order to 
precisely match the upper mass range for SNe II-P. 

(2) The mass range of SNe II-L, albeit narrow, is en- 
tirely consistent with known progenitors of this class men- 
tioned earlier. 

(3) Figure [7] allows SNe Iln to arise from among the 
most massive stars, consistent with their hypothesised LBV 
or pulsational pair instability progenitors. As in the Binary 
#1 and #2 scenarios, it provides for the apparent continuity 
in pre-SN mass loss from SNe II-P to II-L to Iln. The initial 
mass range of SNe Iln progenitors is roughly 23-150 M©, 
commensurate with the known initial mass range of LBVs 
(Smith et al. 2004). We show an alternative version of a 
hybrid scenario in Figure [HI wherein we separate SNe Iln 
and single-star SNe Ic by mass, instead of dividing them 
half-and-half across all single-star masses above 23 M0 . This 
is very similar in principle to the original standard single- 
star hypothesis (Figure [4|, but with SNe lib, lb, and some 
Ic now excluded as binaries. In Figure [H the dividing mass 
between SNe Iln and single-star SNe Ic is ~36 M© . This has 
the advantage that classical LBV eruptions above this mass 
can account for the mass loss to produce SNe Ic, but it has 
the disadvantages that it does not allow SNe Iln to arise 
from the most massive stars, and it does not allow for other 
factors like luminous SNe Iln preferring low metallicity, or 
rapid rotation working across a range of masses. For these 
reasons, we tend to favour Figure [7] over Figure [HI but the 
truth may be somewhere in between. Differentiating between 
these two possibilities is difficult, since we do not yet know 
how to distinguish single-star from binary SNe Ic. 

(4) As in the Binary #2 scenario, SNe lib arise in bina- 
ries, consistent with the progenitor of SN 1993 J (see above). 



The initial mass range of SN lib progenitors in this scenario, 
if they occupy the low-mass end of RLOF binaries, would be 
8.5-12 M0. This is admittedly quite low, and perhaps lower 
than expected for the progenitors of SN 1993J (~15 M0; 
Young et al. 2006) and the SN lib that gave rise to Gas A 
(Krause et al. 2008; Rest et al. 2008), given the strong N 
enrichment in its CSM (Chevalier & Kirshner 1978; Fesen 
& Becker 1991; Chevalier & Oishi 2003). An alternative in- 
terpretation may be that initial rotation rates, metallicity, 
or especially binary separation also play a role here, so that 
some of the SN lib and SN lb progenitors overlap in mass 
range up to 25 Mq depending on these conditions. The wider 
mass range would allow more diversity in the progenitors 
of SNe lib, consistent with the expectations of Chevalier & 
Soderberg (2010). Still, studies thus far have revealed no sur- 
viving companion star for Cas A (Thorstensen et al. 2001; 
Krause et al. 2008), so there may be exceptions where some 
massive single stars produce SNe lib as well. On the other 
hand, we note that Podsiadlowski et al. (1992) expect cases 
where the original secondary star that gains mass in RLOF 
may experience accelerated evolution and explode first, leav- 
ing a widowed SN lib or SN lb progenitor to explode as an 
apparently single stripped-envelope star. Perhaps something 
like this occurred in Cas A. 

(5) The hybrid scenario gives an appealing explanation 
for the tiny observed differences between SNe lib and lb 
(Elmhamdi et al. 2006; Chornock et al. 2010), as in the 
Binary #2 scenario. The initial masses corresponding to 
SNe lb (and SNe Ibc-pec) would then be roughly 12-25 Mq . 
These are massive stars in binaries whose winds can get rid 
of the remaining H, but are not strong enough to fully re- 
move the He envelope, probably because they are underlu- 
minous after RLOF. The SNe lb progenitors likely corre- 
spond to a population of lower luminosity, early-type WN 
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stars that are difficult to detect next to their overluminous 
mass-gainer companions. Perhaps these post-RLOF systems 
would appear as peculiar Be or B[e]-like stars (mainly due to 
their overluminous H-rich companions) in nearby galaxies, 
likely showing signs of asymmetric CSM. 

(6) It retains the quality that SNe Ic will still trace the 
most massive stars, especially those at higher metallicity, 
whether they arise from binaries or single stars. It also gives 
two difi'erent channels for making SNe Ic, perhaps providing 
an avenue for explaining the diversity among SNe Ic (i.e., 
normal vs. broad-lined SNe Ic). This is an important point 
beyond the scope of this paper, but Figure [7] suggests some 
interesting possibilities. Even some broad-lined SNe Ic, how- 
ever, appear to arise from only moderately massive stars, 
based on the ejecta mass estimates and progenitor limits 
(Iwamoto et al. 1994; Mazzah et al. 2002; Sauer et al. 2006; 
Crockett et al. 2007). 

(7) SNe lb, on the other hand, do not trace the highest- 
mass stars or regions of high metallicity quite as well in 
this scenario, since it is probably the lower mass stars or 
lower metallicity stars that fail to drive away their He en- 
velopes. This scenario would predict noticeable differences 
between the environments and progenitors of SNe lb and Ic, 
with SNe Ic tending to trace higher initial mass and higher 
metallicity. There is some empirical support for this (Kelly 
et al. 2008; Anderson & James 2009; Papers I and II), but 
further study should treat SNe lb and Ic separately. 

(8) SNe lb are less common and there are fewer well- 
studied examples compared to SNe Ic, but a recent detailed 
investigation of the SN lb 2007Y revealed a small ejecta 
mass that suggested a low initial mass of only 10-13 Mq 
for the progenitor, and interestingly, deduced a progenitor 
mass-loss rate of only < 10~® Mq (Stritzinger et al. 2009). 
This mass-loss rate derived from radio and X-ray data is 
quite low compared to mass-loss rates of classical WR stars, 
supporting the idea that SNe lb arise from lower-mass stars 
than classical WR stars, and that they have relatively low 
luminosity and weak winds (see also Filippenko 1991). It 
is even possible, for instance, that the wind of the mass- 
gainer companion (e.g., an OB supergiant) will be stronger 
than the wind of the SN lb progenitor star, and that the SN 
blast wave will interact mostly with its companion's wind. 
Whether or not the wind is H-poor is difficult to ascertain 
from radio or X-ray observations, and deriving a progenitor 
mass-loss rate depends also on an assumed wind velocity 
(i.e., it may be significantly lower for a slow B-supergiant 
wind than for a fast WR wind). 

(9) There may be a regime where SNe Ic and SNe Iln 
overlap, coming from the transition between very massive 
single stars that are successful in shedding their H envelopes 
through LBV eruptions (Smith & Owocki 2006) and those 
that cannot. This may depend on initial rotation or metal- 
licity, and we speculate that the transition may be the origin 
of some of the unusual "hybrid" SNe that have been classi- 
fied as Type la/IIn, such as SNe 2002ic, 2005gj, 1997cy, and 
1999E (Turatto et al. 2000; Germany et al. 2000; Hamuy et 
al. 2003; Wood-Vasey et al. 2003; Rigon et al. 2003; Wang 
et al. 2004; Chugai & Yungelson 2004; Kotak et al. 2004; 
Aldering et al. 2006; Benetti et al. 2006; Chugai & Chevalier 
2006; Prieto et al. 2007). Benetti et al. (2006) have argued 
that these may in fact be SNe Ic that appear as SNe Iln be- 
cause of CSM interaction, rather than SNe la; this point is 



speculative and still debated, however. We conjecture that 
unusual SNe Ibn like SN 2006jc (e.g., Pastorello et al. 2007; 
Foley et al. 2007) may fit in a similar transitional category 
of very massive stars. 

All things considered, we favour a hybrid scenario like 
Figure [7] as the basic explanation for the observed fractions 
of various SN types in large galaxies, invoking binary RLOF 
to account for most SNe lib, lb, and some Ic, and yet retain- 
ing single-star mass loss with increasing mass to account for 
SNe II-P, II-L, Iln, and some Ic in the most extreme cases. 
We stress, however, that this is hypothetical, with specific 
binary fractions and other parameters adopted to encapsu- 
late only the broad properties of various SN types. Figure [7] 
adopted a constant fraction of progenitors that go through 
RLOF, whereas this may obviously depend on initial mass, 
and RLOF efficiency may depend on other factors like bi- 
nary separation and metallicity. Thus, the mass divisions 
between various types are meant as a general guide, rather 
than definitive values. This is certainly an oversimplification, 
and there may well be exceptions for individual cases or ex- 
treme conditions. More study is needed, including detailed 
population synthesis models with both binary evolution and 
LBV-like mass loss for massive stars. The binary fraction 
and its variation with initial mass are key parameters, as is 
the behaviour of wind and eruptive mass loss with metallic- 
ity and rotation. However, we hope that keeping a scenario 
such as Figure [7] in mind will be useful to guide intuition for 
mapping SNe to stellar initial masses. 



5 CONSEQUENCES AND FUTURE TESTS 

If SNe lib really result from a different channel than other 
SNe II, a simple comparison of the relative numbers of 
SNe Ibc and SNe II (including SN lib with other SNe II) is 
probably misleading. Such a comparison would make sense 
in the standard single-star scenario where all stars above 
some threshold mass, Mwr, make WR stars and SNe Ibc 
(Figure 31, but we have argued that this simple hypothe- 
sis is contradicted by SN observations. Instead, an analysis 
that retains SN types in-line with the separate binary and 
single-star channels discussed here would be more appro- 
priate. For example, whereas envelope stripping via binary 
RLOF should not necessarily depend on metallicity or initial 
mass (unless the close binary fraction changes with mass), 
the transition SN lib — >■ lb — Ic is caused directly by the 
line-driven wind of the post-RLOF WR-like star (i.e., pro- 
ceeding from a low-luminosity WN with some H, to a normal 
WN, to WC). We should therefore be very interested to see 
how fractions of these subtypes change with metallicity. 

Thus, previous studies that have compared the ratio 
Nibc/Nii, lumping SNe lib together with other SNe II, may 
produce somewhat misleading trends and may inspire erro- 
neous conclusions. In future studies, as larger numbers of all 
types of SNe become available, it will be useful to compare 
relative numbers of individual subtypes (IIb:Ib:Ic) as well as 
the ratio of larger groups that represent different channels 
[e.g., (lib + lh + Ic) / (ILP + H-L + Iln)] with metallic- 
ity and host-galaxy environment. The properties of SN II 
environments would be particularly interesting; the Binary 
#2 or hybrid hypotheses would predict, for example, that 
SNe Iln come from more massive stars and should there- 
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fore trace clusters and H ll regions to a higher degree than 
SNe II-P and II-L. This would not be so noticeable for the 
single-star scenario shown in Figure [H We would not neces- 
sarily expect, however, that SNe Iln would be concentrated 
in galaxy centers, as that may betray a high- met allicity ef- 
fect, which leads instead to SNe Ic for the most massive 
stars. Very massive stars that retain their H envelopes un- 
til shortly before core collapse might instead favour lower 
metallicity, and hence smaller host galaxies. This does in- 
deed seem to be the case, as we point out in Paper II. 

A central hypothesis is that SNe lb trace a population 
of moderately massive stars that have lost their H envelopes 
primarily via binary RLOF. These progenitors are like clas- 
sical WR stars in that they are H deficient, but they differ 
in that they are likely to be underluminous with relatively 
weak winds, and stem from a lower range of initial masses 
of roughly 12-25 Mq. These may not be recognised as WR 
stars because of their weaker winds and less prominent emis- 
sion lines (see also Filippenko 1991). In nearby stellar pop- 
ulations, the SN lb progenitors may be among the group 
of underluminous early WN stars, or they may reside in bi- 
nary systems where they are hard to detect next to their 
overluminous mass-gainer companions. We have speculated 
that these post-RLOF systems of moderate mass may ap- 
pear as Be or B[e]-like stars, perhaps with asymmetric CSM. 
Other potential SN lb or Ic progenitor systems are famous 
WR-hOB systems like V444 Cygni, 7^ Vel, or RY Set. 

Although SNe lb are relatively rare, it will be important 
to distinguish SNe lb from Ic in future analyses, and to clar- 
ify their different properties as well as any range in param- 
eter space where they may overlap. It will be especially im- 
portant to further clarify the residual He surface mass that 
separates SNe lb from Ic; there may obviously be examples 
of a smooth transition in He mass between them. If SNe lb 
arise from RLOF in binary systems, then the mass-loss rates 
of the progenitor stars derived from radio and X-ray obser- 
vations may be tricky to interpret. For example, we noted 
that the wind of the overluminous mass-gainer companion 
may be stronger than the SN progenitor star itself, and so 
interaction between the SN blast wave and the companion's 
wind might dominate the observed radio and X-ray emis- 
sion. Without a radiative shock to produce strong Balmer 
lines, it would be difficult to determine whether the wind is 
deficient in hydrogenQ 

Lastly, it would be interesting to further investigate 
differences among environments and progenitors of SNe Ic, 
since this class alone makes up 15% of all CCSNe (a substan- 
tial fraction of the most massive stars), and may have mul- 
tiple progenitor channels. Do the broad-lined SNe Ic arise 
preferentially from one channel? This is a key question in re- 
gard to the progenitors of long-duration gamma-ray bursts. 



^ In the special case of SN 2006jc, the dense CSM produced 
strong He I lines, so one can infer that the progenitor star suffered 
a precursor eruption and that the CSM was not from a compan- 
ion (Pastorello ct al. 2007; Foley et al. 2007; Smith et al. 2008a); 
other cases of SNe lb are less clear. 



6 CONCLUSIONS 

We have studied the observed fractions of different SN types 
from LOSS, and considered the implications for massive 
star evolution. Assuming a Salpeter IMF, we have exam- 
ined what ranges of initial mass are needed to account for 
the observed fractions of SNe II-P, II-L, Iln, lib, lb, and Ic 
under various assumptions about the roles of stellar winds 
and close binary RLOF in stellar evolution. We briefiy list 
the main conclusions here, which apply to stellar evolution 
in relatively large galaxies. 

(1) A major finding is that the high observed fraction of 
SNe Ibc cannot be reconciled with predictions of single-star 
evolution, where a star's own wind dominates the removal 
of its H envelope. The initial-mass range corresponding to 
the observed population of classical WR stars can only ac- 
count for about half of the observed SNe Ibc, so classical 
WR stars are not the progenitors of a significant fraction of 
SNe Ibc. Similarly, the initial mass above which single stars 
are expected to shed their own envelopes provides a vastly 
insufficient fraction of stripped-envelope progenitors, even 
with the overly generous mass-loss rates adopted in most 
stellar evolution models. 

(2) Instead, we find it likely that RLOF in binary sys- 
tems is responsible for the stripped-envelope progenitors of 
most SNe lib and lb, and probably a large fraction of SNe Ic 
as well. If these are distributed over the full range of masses, 
then SNe lib and lb probably arise from lower initial masses 
of 8.5-25 Mq , and SNe Ic arise from more massive stars with 
stronger winds. 

(3) Even if binary RLOF dominates the removal of the 
H envelope, the further removal of the He layer depends on 
metallicity-dependent line-driven winds of the WR star, so 
SNe Ic are still expected to favour more luminous stars and 
higher metallicity environments. 

(4) If the progenitors of SNe lb and lib are not classical 
WR stars because their initial masses are too low, then what 
kind of stars are the progenitors? We conjecture that they 
are probably underluminous H-poor stars with weak winds 
that would not necessarily be recognised as WR stars with 
prominent emission-line spectra. They may be easily hidden 
by their overluminous mass-gainer companions, which may 
in some cases appear as B[e] supergiants or related stars 
with asymmetric CSM. If so, one must be cautious when 
interpreting signatures of the CSM interaction in SNe lb, as 
the emission may in some cases be dominated by SN shock 
interaction with a companion star's wind. 

(5) If binary RLOF is important in producing stripped- 
envelope progenitors that are a substantial fraction (1/4 to 
1/3) of all SN progenitors, then it would be a mistake to use 
statistics of SN types or WR/O-star ratios to guide models 
for single-star evolution. 

(6) After shifting most stripped-envelope progenitors to 
the binary RLOF channel, the progenitors of the remain- 
ing SN types (H-rich single stars, wide binaries, and possi- 
bly mass-gainers in RLOF binaries) must be redistributed 
across the full range of initial masses. In our favoured sce- 
nario (Figure [T]), SNe II-P correspond to initial masses of 
roughly 8.5-18 Mq, SNe H-L to 18-23 Mq, and SNe Iln to 
23-150 Mq . This produces good agreement with mass ranges 
inferred from progenitor studies of SNe II-P, II-L, and Iln. 
In particular, this allows some SNe Iln to arise from among 
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the most massive stars, as suggested by some very luminous 
SNe Iln. Most stellar evolution models fail to account for 
very massive stars reaching core collapse without shedding 
their H envelope, but this is an expected outcome of the 
lower mass- loss rates now dictated by observations. We also 
find it likely that some fraction of the most massive sin- 
gle stars shed their H envelopes to produce WR stars and 
SNe Ic, probably due to high metallicity. This allows for the 
possibility that SNe Iln favour low- metallicity environments. 

(7) We briefly consider the possibility that some massive 
stars collapse directly to BHs without a visible SN display. 
We can rule out this option for the scenario of standard 
single-star evolution, because it would make all the prob- 
lems we note with Figure|4]worse. We find no empirical sup- 
port for the argument that the "RSG problem" may imply 
direct SN-less BH formation, because this problem largely 
goes away with more reliable SN subtype fractions and with 
the realization that some RSG stars evolve to other types 
of progenitors before exploding. Though we cannot rule out 
the possibility that some massive single stars within a par- 
ticular mass range suffer quiet collapse to a BH in a binary- 
dominated scenario, quiet BIf collapse is not required to 
explain the observed relative fractions of CCSNe. 

ACKNOWLEDGMENTS 

We are grateful to the many students, postdocs, and other collabora- 
tors who have contributed to the Katzman Automatic Imaging Telescope 
and the Lick Observatory Supernova Search over the past two decades. We 
thank the Lick Observatory staff for their assistance with the operation of 
KAIT. A.V.F.'s group has been supported by many grants from the US 
National Science Foundation (NSF; most recently AST-0607485 and AST- 
0908886), the TABASCO Foundation, US Department of Energy SciDAC 
grant DE-FC02-06ER41453, and US Department of Energy grant DE-FG02- 
08ER41653. KAIT and its ongoing operation were made possible by dona- 
tions from Sun Microsystems, Inc., the Hewlett-Packard Company, Auto- 
Scope Corporation, Lick Observatory, the NSF, the University of California, 
the Sylvia &c Jim Katzman Foundation, and the TABASCO Foundation. 



REFERENCES 

Aldering, G., Humphreys, R.M., & Richmond, M. 1994, AJ, 107, 
662 

Aldering, G., et al. 2006, ApJ, 650, 510 
Anderson, J. P., & James, P.A. 2009, MNRAS, 339, 559 
Arcavi, I., et al. 2010, ApJ, submitted l larXiv:1004.0615l 
Arnett, W.D. 1987, ApJ, 319, 136 

Arnett, W.D., Bahcall, J.N., Kirshner, R.P., & Woosley, S.P. 

1989, ARA&A, 27, 629 
Barth, A.J., Van Dyk, S.D., Filippenko, A.V., Leibundgut, B., & 

Richmond, MW. 1996, AJ, 111, 2047 
Bastian, N., Covey, K.R., & Meyer, M.R. 2010, submitted 

jarXiy;l 001. 2965 ) 
Benetti, S., Cappellaro, E., Turatto, M., Taubenberger, S., 

Haratyunyan, A., & Valenti, S. 2006, ApJ, 654, L12 
Boissier, S., & Prantzos, N. 2009, A&A, 503, 137 
Bouret, J.C., Lanz, T., & Hillier, D.J. 2005, A&A, 438, 301 
Branch, D., Nomoto, K., & Filippenko, A.V. 1991, Comm. As- 

troph., 15, 221 

Brown, T.M., Heap, S.R., Hubeny, I., Lanz, T., & Lindler, D. 

2002, ApJ, 579, L75 
Cappellaro, E., Turatto, M., Tsvetkov, D.Y., Bartunov, O.S., Pol- 
las, C, Evans, R., & Hamuy, M. 1997, A&A, 322, 431 
Chevalier, R.A., & Kirshner, R.P. 1978, ApJ, 219, 931 
Chevalier, R.A., & Oishi, J. 2003, ApJ, 593, L23 
Chevalier, R.A., & Soderberg, A. 2010, ApJ, 711, L40 
Chornock, R., et al. 2010, submitted llarXiv:1001. 277511 
Chugai, N.N. 1991, MNRAS, 250, 513 



Chugai, N.N., & Chevalier, R.A. 2006, ApJ, 641, 1051 

Chugai, N.N., & Yungolson, L.R. 2004, Astron. Lett., 30, 83 

Chugai, N.N., et al. 2004, MNRAS, 352, 1213 

Conti, P.S. 1976, Societe Royale des Sciences de Liege, 9, 193 

Crockett, R.M., et al. 2007, MNRAS, 381, 835 

Crowther, P.A. 2007, ARA&A, 45, 177 

Crowther, P.A., & Hadficld, L. 2006 

Crowther, P.A., et al. 2006 

Davidson, K., et al. 1982, ApJ, 253, 696 

de Jager, C, Nieuwenhuijzen, H., & van der Hught, K.A. 1998, 

A&AS, 72, 259 
Doggett, J.B., & Branch, D. 1985, AJ, 90, 2303 
Eldridge, J.J., & Tout, C.A. 2004, MNRAS, 353, 87 
Eldridge, J. J., Izzard, R.G., & Tout, C.A. 2008, MNRAS, 384, 

1109 

Elias-Rosa, N., Van Dyk, S.D., Li, W., Morrell, N., Gonzalez, S., 
Hamuy, M., Filippenko, A.V., Cuillandre, J.C., Foley, R., & 
Smith, N. 2010a, ApJ, 706, 1174 

Elias-Rosa, N., et al. 2010b, in press llarXiv:0912.2880ll 

Elias-Rosa, N., et al. 2010c, submitted 

Elmhamdi, A., Danziger, LJ., Branch, D., Leibundgut, B, Baron, 

E., & Kirshner, R.P. 2006, A&A, 450, 305 
Falk, S.W., & Arnett, W.D. 1977, ApJS, 33, 515 
Fesen, R.A, & Becker, R.H. 1991, ApJ, 371, 621 
Figer, D.F. 2005, Nature, 434, 192 
Filippenko, A. V. 1991, lAU Symposium 143, 529 
Filippenko, A. V. 1997, ARAA, 35, 309 

Filippenko, A.V., Matheson, T., & Barth, A.J. 1994, AJ, 108, 
2220 

Filippenko, A.V., Li, W.D., Treffers, R.R., & Modjaz, M. 2001, in 
Small- Telescope Astronomy on Global Scales, ed. W.P. Chen, 

C. Lemme, & B. Paczynski (San Francisco; ASP), 121 
Foley, R.J., Smith, N., Ganeshalingam, M., et al. 2007, ApJ, 657, 

L105 

Eraser, M., et al. 2010, in press larXiv:0912.2071 ll 
Fryer, C.L. 1999, ApJ, 522, 413 

Fryer, C.L., Burrows, A., & Bez, W. 1998, ApJ, 496, 333 
Fryer, C.L., Woosley, S.P., & Hartmann, D.H. 1999, ApJ, 526, 
152 

PuUerton, A.W., Massa, D.L., & Prinja, R.K. 2006, ApJ, 637, 
1025 

Gal- Yam, A., & Leonard, D.C. 2009, Nature, 458, 865 
Garnett, D.R. 2002, ApJ, 581, 1019 

Germany, L.M., Reiss, D.J., Sadler, E.M., Schmidt, B.P., & 

Stubbs, C.W. 2000, ApJ, 533, 320 
Hamuy, M., et al. 2003, Nature, 424, 651 

Heger, A., Fryer, C.L., Woosley, S.E., Langer, N., & Hartmann, 

D. H. 2003, ApJ, 591, 288 

Hendry, M.A., et al. 2006, MNRAS, 369, 1303 
Humphreys, R.M., & McElroy, D.B. 1984, ApJ, 284, 565 
Humphreys, R.M., Nichols, M., & Massey, P. 1985, AJ, 90, 101 
Iwamoto, K., Nomoto, K., Hofiich, P., Yamaoka, H., Kumagai, S., 

& Shigeyama, T. 1994, ApJ, 437, L1151 
Izotov, Y.I., Foltz, C.B., Green, R.F., Guseva, N.G., & Thuan, 

T.X. 1997, ApJ, 487, L37 
Kelly, P.L., Kirshner, R.P., & Pahre, M. 2008, ApJ, 687, 1201 
Kitaura, F.S., Janka, H.T., & Hillebrandt, W. 2006, A&A, 450, 

345 

Kobulnicky, H.A., & Fryer, C.L. 2007, ApJ, 670, 747 
Kochanek, C.S., Beacom, J.F., Kistler, M.D., Prieto, J.L., & 

Stanek, K.Z. 2008, ApJ, 684, 1336 
Kotak, R., & Vink, J.S. 2006, A&A, 460, L5 
Kotak, R., et al. 2004, MNRAS, 354, L13 
Krause, O., et al. 2008, Science, 320, 1195 
Kudritzki, R.P., & Puis, J. 2000, ARAA, 38, 613 
Lamers, H.J.G.L.M., & Cassinelli, J. P. 1999, Introduction to Stel- 
lar Winds (Cambridge: Cambridge University Press) 
Leaman, J., et al. 2010, MNRAS, in prep. (Paper I) 



18 N. Smith et al. 



Leloudas, G., SoUcrman, J., Levan, A. J., Fynbo, J.P.U., Malesani, 

D., & Maund, J.R. 2010, A&A, in press llarXiv:1002.3164 ) 
Leonard, D.C., et al. 2002, AJ, 124, 2490 

Leonard, D.C., Kanbur, S.M., Ngeow, C.C., & Tanvir, N.R. 2003, 
ApJ, 594, 247 

Li, W., Van Dyk, S.D., Filippenko, A.V., & Cuillandre, J.C. 2005, 
PASP, 117, 121 

Li, W., Van Dyk, S.D., Filippenko, A.V., Cuillandre, J.C, Jha, 
S., Bloom, J.S., Reiss, A.G., & Livio, M. 2006, ApJ, 641, 1060 
Li, W., Wang, X., Van Dyk, S.D., Cuillandre, J.C, Foley, R., & 

Filippenko, A.V. 2007, ApJ, 661, 1013 
Li, W., et al. 2010a, MNRAS, in prep. (Paper II) 
Li, W., et al. 2010b, MNRAS, in prep. (Paper III) 
Litvinova, I.Y., & Nadyozhin, D.K. 1983, Ap. Sp. Sci., 89, 89 
Massey, P. 2003, ARAA, 41, 15 

Massey, P., Johnson, K.E., k, DeCioia-Eastwood, K. 1985, ApJ, 
454, 151 

Massey, P., DeCioia-Eastwood, K., & Waterhouse, E. 2001, A J, 
121, 1050 

Maund, J.R., & Smartt, S.J. 2005, MNRAS, 360, 288 

Maund, J.R., & Smartt, S.J. 2009, Science, 324, 486 

Maund, J.R., Smartt, S.J., Kudritzki, R.P., Podsiadlowski, P., & 

Cilmore, C.F. 2004, Nature, 427, 129 
Maund, J.R., Smartt, S.J., & Danziger, I.J. 2005, MNRAS, 364, 

L33 

Maund, J.R., Wheeler, J.C, Ferdinando, P., Wang, L., Baade, 

D., & Hoflich, P.A. 2007, ApJ, 671, 1944 
Mazzah, P.A., et al. 2002, ApJ, 572, L61 
Meszaros, Sz., Avrett E.H., & Dupree, A. 2009, AJ, 138, 615 
Meynet, C, Maeder, A., Schaller, C, Shaerer, D., & Charbonnel, 

C. 1994, A&A Supp., 103, 97 
Meynet, C, Mowlvi, N., & Maeder, A. 2008, preprint 

(astro-ph/0611261 ) 
Miller, A. A., et al. 2009, ApJ, 690, 1303 
Modjaz, M., et al. 2008, AJ, 135, 1136 

Nelemans, C, Voss, R., Nielsen, M.T.B., & Roelofs, C. 2010, 

MNRAS, in press (arXiv: 10021657) 
Nieuwenhuijzen, H., & de Jager, C. 1990, A&A, 231, 134 
Nomoto, K. 1984, ApJ, 277, 791 
Nomoto, K., et al. 1993, Nature, 364, 507 
Paczyhski, B. 1967, Acta Astron., 17, 355 
Pastorello, A., et al. 2007, Nature, 449, 1 

Podsiadlowski, P., Joss, P.C, & Hsu, J.J.L. 1992, ApJ, 391, 246 
Podsiadlowski, P., Hsu, J.J.L., Joss, P.C, & Ross, R.R. 1993, 

Nature, 364, 509 
Prantzos, N., & Boissier, S. 2003, A&A, 406, 259 
Prieto, J.L., et al. 2007, submitted ( arXiv:0706.4088 ) 
Prieto, J.L., Stanek, K.Z., & Beacom, J.F. 2008, ApJ, 673, 999 
Puis, J., et al. 2006, A&A, 454, 625 
Pumo, M.L., et al. 2009, submitted (arXiv:0910.0640j) 
Quimby, R., et al. 2009, submitted (arXiv:0910:0059) 
Reimers, D. 1977, A&A, 61, 217 
Rest, A., et al. 2008, ApJ, 681, L81 
Rigon, L., et al. 2003, MNRAS, 340, 191 

Ryder, S.D., Murrowood, C.E., & Stathakis, R.A. 2006, MNRAS, 
369, L32 

Ryder, S.D., Sadler, E.M., Subrahmanyan, R., Weiler, K.W., 

Panagia, N., & Stockdale, C. 2004, MNRAS, 349, 1093 
Salpeter, E.E. 1955, ApJ, 121, 161 
Sauer, et al. 2006, MNRAS, 369, 1939 

Schinzel, F.K., Taylor, G.B., Stockdale, C.J., Granot, J., & 

Ramirez-Ruiz, E. 2009, ApJ, 691, 1380 
Silverman, J., et al. 2009, PASP, 121, 689 
Smartt, S.J. 2009, ARA&A, 47, 63 

Smartt, S.J., Eldridge, J.J., Crockett, R.M., & Maund, J.R. 2009, 

MNRAS, 395, 1409 
Smartt, S.J., Cilmore, G.F., Trentham, N., Tout, C.A., & Frayn, 

CM. 2001, ApJ, 556, L29 



Smartt, S.J., Cilmore, C.F., Tout, C.A., & Hodgkin, S.T. 2002, 

ApJ, 565, 1089 
Smartt, S.J., et al. 2003, MNRAS, 343, 735 
Smartt, S.J., et al. 2004, Science, 303, 499 
Smith, N. 2007, AJ, 133, 1034 

Smith, N. 2008, lAU Symposium 250: Massive Stars as Cosmic 
Engines, eds. F. Bresolin, P.A. Crowther, & J. Puis (Cam- 
bridge: Cambridge Univ. Press), 193 

Smith, N., Chornock, R., Silverman, J.M., FiUppcnko, A.V., & 
Foley, R.J. 2010, ApJ, 709, 856 

Smith, N., & Conti, P. S. 2008, ApJ, 679, 1467 

Smith, N., Foley, R.J., & Filippenko, A.V. 2008a, ApJ, 680, 568 

Smith, N., Hinkle, K.H., & Ryde, N. 2009a, AJ, 137, 3558 

Smith, N., & McCray, R. 2007, ApJ, 671, L17 

Smith, N., & Owocki, S. P. 2006, ApJ, 645, L45 

Smith, N., Vink, J., & de Koter, A. 2004, ApJ, 615, 475 

Smith, N., et al. 2007, ApJ, 666, 1116 

Smith, N., et al. 2008b, ApJ, 686, 467 

Smith, N., et al. 2009b, ApJ, 695, 1334 

Smith, N., et al. 2010, in preparation 

Sramek, R.A., & Weiler, K.W. 1990, in Supernovae, ed. A.C. 

Petschek (New York: Springer), 76 
Stanek, K.Z., et al. 2006, Acta Astronomica, 56, 333 
Stritzinger, M., et al. 2009, ApJ, 696, 713 

Thompson, T.A., Prieto, J.L., Stanek, K.Z., et al. 2009, ApJ, 705, 
1364 

Thorstensen, J.R., Fesen, R.A., & van den Bergh, S. 2001, AJ, 
122, 297 

Trundle, C, Kotak, R., Vink, J.S., & Meikle, W.P.S. 2008, A&A, 
483, L47 

Turatto, M., et al. 2000, ApJ, 534, L57 

Vanbeveren, D., Van Bever, J., & Belkus, H. 2007, ApJ, 662, L107 
Van den Bergh, S., McClure, R.D., & Evans, R. 1987, ApJ, 323, 
44 

Van Dyk, S.D., Peng, C.Y., Varth, A.J., & FiUppcnko, A.V. 1999, 
AJ, 118, 2331 

Van Dyk, S.D. 2010, in Eta Carinae and the Supernova Impostors, 
ed. R.M. Humphreys & K. Davidson (New York: Springer), in 
press 

Van Dyk, S.D., Li, W., & Filippenko, A.V. 2003a, PASP, 115, 1 
Van Dyk, S.D., Li, W., & Filippenko, A.V. 2003b, PASP, 115, 448 
Van Dyk, S.D., Li, W., & Filippenko, A.V. 2003c, PASP, 115, 
1289 

Van Dyk, S.D., et al. 2002, PASP, 114, 1322 
Van Dyk, S.D. et al. 2009, AAS Bull., 214, 604.02 
Wanajo, S.E., et al. 2009, ApJ, 695, 208 
Wang, L., et al. 2004, ApJ, 604, L53 
Wellstein, S., & Langer, N. 1999, A&A, 350, 148 
Wheeler, J.C, & Swartz, D.A. 1993, Sp. Sci. Rev., 66, 425 
Wood-Vasey, W.M., Wang, L., & Aldering, C. 2004, ApJ, 616, 
339 

Woosley, S.E., Blinnikov, S., & Heger, A. 2007, Nature, 450, 390 
Woosley, S.E., & Bloom, J.S. 2006, ARA&A, 44, 507 
Woosley, S.E., Eastman, R.G., Weaver, T.A., & Pinto, P.A. 1994, 
ApJ, 429, 300 

Woosley, S.E., Heger, A., & Weaver, T.A. 2002, Rev. Mod. Phys., 

74, 1015 

Yoon, S.C, & Cantiello, M. 2010. larXiv: 1005.49251 
Young, D.R., et al. 2006, ApJ, 640, 891 
Young, D.R., et al. 2008, A&A, 489, 359 



